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ABSTRACT 

We investigate the properties of the galaxies selected from the deepest 850-^m survey 
undertaken to date with SCUBA-2 on the JCMT as part of the SCUBA-2 Cosmology 
Legacy Survey. A total of 106 sources (>5 (t) were uncovered at 850 ^m from an area 
of ~ 150 arcmin^ in the centre of the COSMOS/UltraVISTA/CANDELS field, imaged 
to a typical depth of CTgso — 0.25 mJy. We utilise the available multi-frequency data to 
identify galaxy counterparts for 80 of these sources (75%), and to establish the com¬ 
plete redshift distribution for this sample, yielding z = 2.38 ± 0.09. We have also been 
able to determine the stellar masses of the majority of the galaxy identifications, en¬ 
abling us to explore their location on the star-formation rate : stellar-mass (SFR:M*) 
plane. Crucially, our new deep 850-/rm selected sample reaches flux densities equiv¬ 
alent to SFR ~ lOOM 0 yr“^, enabling us to confirm that sub-mm galaxies form the 
high-mass end of the ‘main sequence’ (MS) of star-forming galaxies at z > 1.5 (with 
a mean specific SFR of sSFR = 2.25 ± 0.19 Gyr“^ at z ~ 2.5). Our results are consis¬ 
tent with no significant flattening of the MS towards high masses at these redshifts. 
However, our results add to the growing evidence that average sSFR rises only slowly 
at high redshift, resulting in logj^psSFR being an apparently simple linear function of 
the age of the Universe. 

Key words: galaxies: high-redshift, evolution, starburst - cosmology: observations - 
submillimetre: galaxies 


1 INTRODUCTION 

It is now well known that approximately half of the starlight 
in the Universe is re-processed by cosmic dust and re¬ 
emitted at far-infrared wavelengths (Dole et al. 2006). How¬ 
ever, due to a combination of the inescapable physics of 
diffraction, the molecular content of our atmosphere, and 
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the technical difficulties of sensitive high-background imag¬ 
ing, it has proved difficult to connect the UV/optical and 
far-infrared/sub-mm views of the Universe into a consis¬ 
tent and complete picture of galaxy formation/evolution. 
Thus, while the advent of SCUBA on the 15-m James Clerk 
Maxwell Telescope (JCMT) in the late 1990s (Holland et al. 
1999) enabled the first discovery of distant dusty galaxies 
with star-formation rates SFR ~ 1000 M 0 yr“^ (e.g. Small 
et al. 1997; Hughes et al. 1998; Barger et al. 1998; Bales 
et al. 1999), such objects initially seemed too extreme and 
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unusual to be easily related to the more numerous, ‘normal’ 
star-forming galaxies being uncovered at UV/optical wave¬ 
lengths at comparable redshifts (« ~ 2 — 4) by Keck (e.g. 
Steidel et al. 1996) and the Hubble Space Telescope (HST) 
(e.g. Madau et al. 1996). In recent years the study of rest- 
frame UV-selected galaxies has been extended out beyond 
2 ~ 10 (see Dunlop 2013 for a review, and Coe et al. 2013; 
Ellis et al. 2013; McLure et al. 2013; Bouwens et al. 2015; 
Bowler et al. 2012, 2014; Oesch et al. 2014; Finkelstein et al. 
2015; Ishigaki et al. 2015; McLeod et al. 2015), while a num¬ 
ber of sub-mm selected galaxies have now been confirmed at 
z > 4 (Capak et al. 2008; Coppin et al. 2009; Daddi et al. 
2009a,b; Knudsen et al. 2010; Riechers et al. 2010; Cox et 
al. 2011; Combes et al. 2012; Weiss et al. 2013) with the 
current redshift record holder at z = 6.34 (Riechers et al. 
2013). However, while such progress is exciting, at present 
there is still relatively little meaningful intersection between 
these UV/optical and far-infrared/sub-mm studies of the 
high-redshift Universe (although see Walter et al. 2012). 

At more moderate redshifts, however, recent years have 
seen increasingly successful efforts to bridge the gap between 
the unobscured and dust-enshrouded views of the evolving 
galaxy population. Of particular importance in this endeav¬ 
our has been the power of deep 24 pm imaging with the 
MIPS instrument on board Spitzer, which has proved ca¬ 
pable of providing a useful estimate of the dust-obscured 
star-formation activity in a significant fraction of optically- 
selected galaxies out to z ~ 1.5 — 2 (e.g. Caputi et al. 2006; 
Elbaz et al. 2010). Indeed, MIPS imaging of the GOODS sur¬ 
vey fields played a key role in establishing what has proved 
to be a fruitful framework for the study of galaxy evolution, 
namely the existence of a so-called “main sequence” (MS) for 
star-forming galaxies, in which star-formation rate is found 
to be roughly proportional to stellar mass (SFR oc M,; 
Noeske et al. 2007; Daddi et al. 2007; Renzini & Peng 2015), 
with a normalisation that rises with increasing redshift (e.g. 
Santini et al. 2009; Oliver et al. 2010; Elbaz et al. 2011; 
Karim et al 2011; Rodighiero et al. 2011, 2014; Tasca et al. 
2015; Salmon et al. 2015; Schreiber et al. 2015; Johnston et 
al. 2015). 

Interest in the MS of star-forming galaxies has contin¬ 
ued to grow (see Speagle et al. 2014 for a useful and com¬ 
prehensive overview), not least because of the difficulty en¬ 
countered by most current models of galaxy formation in re¬ 
producing its apparently rapid evolution between z ~ 0 and 
« ~ 2 (e.g. Mitchell et al. 2014). However, it has, until now, 
proved very difficult to extend the robust study of the MS 
beyond z cz 2 and to the highest masses (e.g. Steinhardt et 
al. 2014; Salmon et al. 2015; Leja et al. 2015). This is because 
an increasing fraction of star formation is enshrouded in dust 
in high-mass galaxies, and Spitzer MIPS and Herschel be¬ 
come increasingly ineffective in the study of dust-enshrouded 
SF with increasing redshift (due to a mix of wavelength 
and resolution limitations), as the far-infrared emission from 
dust is redshifted into the sub-mm/mm regime. 

A complete picture of star-formation in more massive 
galaxies at high-redshift can therefore only be achieved with 
ground-based sub-mm/mm observations, which provide im¬ 
age quality at sub-mm wavelengths that is vastly superior 
to what can currently be achieved from space. The chal¬ 
lenge, then, is to connect the population of dusty, rapidly 
star-forming high-redshift galaxies revealed by ground-based 


sub-mm/mm surveys to the population of more moderate 
star-forming galaxies now being revealed by optical/near- 
infrared observations out to the highest redshifts. On a 
source-by-source basis this can now be achieved by targeted 
follow-up of known optical/infrared-selected galaxies with 
ALMA (e.g. Ono et al. 2014). However, this will inevitably 
produce a biased perspective which can only be re-balanced 
by also continuing to undertake ever deeper and wider sub- 
mm/mm surveys capable of detecting highly-obscured ob¬ 
jects (again, potentially, for ALMA follow-up; Karim et al. 
2013; Hodge et al. 2013), and thus completing our census of 
star-forming galaxies in the young Universe. 

This is one of the primary science drivers for the 
SCUBA-2 Cosmology Legacy Survey (S2CLS). The S2CLS 
is advancing the field in two directions. First, building on 
previous efforts with SCUBA (e.g. Scott et al. 2002; Scott, 
Dunlop & Serjeant 2006; Coppin et al. 2006), MAMBO 
(e.g. Bertoldi et al. 2000; Greve et al. 2004), LABOCA 
(e.g. Weiss et al. 2009) and AzTEC (e.g. Austermann et 
al. 2010; Scott et al. 2012), the S2CLS is using the im¬ 
proved mapping capabilities of SCUBA-2 (Holland et al. 
2013) to extend surveys for bright (Sgso > 5 mjy) sub- 
mm sources to areas of several square degrees, yielding 
large statistical samples of such sources (> 1000). Second, 
the S2CLS is exploiting the very dryest (Grade-1) condi¬ 
tions at the JCMT on Manna Kea, Hawaii to obtain very 
deep 450 pm imaging of small areas of sky centred on the 
HST CANDELS fields (Grogin et al. 2011), which provide 
the very best multi-wavelength supporting data to facili¬ 
tate galaxy counterpart identification and study. The first 
such deep 450 pm image has been completed in the centre 
of the COSMOS-CANDELS/UltraVISTA field, with the re¬ 
sults reported by Geach et al. (2013) and Roseboom et al. 
(2013). Here we utilise the ultra-deep 850/rm image of the 
same region, which was automatically obtained in parallel 
with the 450 pm imaging. While the dryest weather is more 
essential for the shorter-wavelength imaging at the JCMT, 
such excellent conditions (and long integrations) inevitably 
also benefit the parallel 850/rm imaging. Consequently, the 
850 pm data studied here constitute the deepest ever 850 pm 
survey ever undertaken over an area ~ 150arcmin^. 

The depth of the new S2CLS 850 pm imaging is typi¬ 
cally (T 850 — 0.25 mJy. This is important because it means 
that galaxies detected near the limit of this survey have 
SFR ~ lOOM 0 yr“^, which is much more comparable to 
the highest SFR values derived from UV/optical/near- 
infrared studies than the typical SFR sensitivity achieved 
with previous single-dish sub-mm/mm imaging (i.e. SFR ~ 
lOOOM 0 yr“^ as a result of crgso — 2mJy). Ultimately, of 
course, ALMA will provide even deeper sub-mm surveys 
with the resolution required to overcome the confusion limit 
of the single-dish surveys. However, because of its modest 
field of view (~ 20 arcsec at 850 pm) it is observationally 
expensive to survey large areas of blank sky with ALMA, 
and contiguous mosaic surveys are hard to justify at depths 
where the source surface density is significantly less than one 
per pointing. Thus, at the intermediate depths probed here, 
the S2CLS continues to occupy a unique and powerful niche 
in the search for dust-enshrouded star-forming galaxies. 

The fact that previous sub-mm/mm surveys were only 
generally capable of detecting very extreme objects has un¬ 
doubtedly contributed to some of the confusion/controversy 
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Figure 1. The SCUBA-2 850 map of the central sub-region of the COSMOS/UltraVISTA field. All 106 sources extracted with 
SNR > 5 (t are highlighted with red circles, and marked with the ID number by which the sources are tabulated in Table [AT] As explained 
in Section |2.2| two optical-infrared catalogues were utilised in this work. The catalogue with the HST CANDELS and deconfused IRAC 
SEDS data, which contains sources extracted from the map covering the CANDELS area (enclosed by the two blue vertical lines), and 
the catalogue without the HST data (outside the blue lines) which was used only for 20 sources. 


over the nature of galaxies selected at sub-mm/mm wave¬ 
lengths; while Michalowski et al. (2012b) and Roseboom et 
al. (2013) have presented evidence that sub-mm selected 
galaxies lie on the high mass end of the MS at 2 = 2 —3, oth¬ 
ers have continued to argue that, like many local ULIRGs, 
they are extreme pathological objects driven by recent ma¬ 
jor mergers (e.g. Mainline et al. 2011). Some of this debate 
reflects disagreements over the stellar masses of the objects 
rather than their star-formation rates (e.g. Michalowski et 
al. 2014). Nevertheless, the fact that even high-mass galax¬ 
ies on the MS lay right at the detection limits of previ¬ 
ous sub-mm surveys inevitably resulted in many sub-mm 
selected objects apparently lying above the MS, fueling ar¬ 
guments about whether they were indeed significant outliers. 


or whether we have simply been uncovering the positive tail 
in SFR around the MS (see Roseboom et al. 2013). 

The much deeper 850 fim survey studied here is capa¬ 
ble of settling this issue, provided of course we can overcome 
the now customary challenge of identifying the galaxy coun¬ 
terparts of most of the sub-mm sources, and determining 
their redshifts, SFRs and stellar masses (M*) (e.g. Ivison et 
al. 2007; Dunlop et al. 2010; Biggs et al. 2011; Wardlow et 
al. 2011; Michalowski et al. 2012a; Koprowski et al. 2014). 
However, in this effort, we are also aided by the depth of the 
SCUBA-2 data, and by the additional positional information 
provided by the (unusual) availability of 450 fim detections 
(with FWHM ~ Sarcsec) for 50% of the sample. We also 
benefit hugely from the unparalleled multi-frequency sup- 
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porting data available in the CANDELS fields, provided by 
HST, Subaru, CFHT, Vista, Spitzer, Herschel and the VLA. 

The paper is structured as follows. In Section 2 we 
present the SCUBA-2 and other multi-wavelength data uti¬ 
lized in this work. Then, in Section 3 we explain how op¬ 
tical/infrared galaxy counterparts were established for the 
SCUBA-2 sources, and summarize the resulting identifica¬ 
tion statistics. Next, in Section 4 we explain the calcula¬ 
tion of the photometric redshifts, both from the optical- 
infrared data for the galaxy identifications, and from the 
long-wavelength data for the unidentified or spuriously iden¬ 
tified sources. The resulting redshift distribution for the 
complete 106-source S2CLS sample is presented here, and 
compared with the redshift distributions derived from other 
recent sub-mm/mm surveys. In Section 5 we move on to 
derive and discuss the physical properties of the sources 
(such as dust temperature, bolometric luminosity, SFR, stel¬ 
lar mass), culminating in the calculation of specific SFR, 
and the exploration of the star-forming MS. Our conclusions 
are summarized in Section 6 . All magnitudes are quoted 
in the AB system (Oke 1974; Oke & Gunn 1983) and all 
cosmological calculations assume Qm ~ 0.3, Ha = 0.7 and 
Ho = 70kms“^Mpc“^. 
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Figure 2. The distribution of the separation between 850 
and 450 /rm positions, where the maximum search radius was set 
to 25 arcsec. Based on this distribution, it was decided that a 
450 /rm source and an 850 /rm source lying within 6 arcsec of each 
other correspond to the same galaxy (for such pairs, the mean 
separation is 2.7 ± 0.2 arcsec). As detailed in Table [aT] there are 
53 850/rm sources with 450 ;rm counterparts (50%). 


2 DATA 

2.1 SCUBA-2 imaging &: source extraction 

We used the deep 850 fim and 450 p.m S2CLS imaging of the 
central ~ 150arcmin^ of the COSMOS/UltraVISTA field, 
coincident with the Spitzer SEDS (Ashby et al. 2013) and 
HST CANDELS (Grogin et al. 2011) imaging. The obser¬ 
vations were taken with SCUBA-2 mounted on the JCMT 
between October 2011 and March 2013, reaching depths of 
(Tsso — 0.25 mjy and (T 450 — 1.5 mjy (Geach et al. 2013, 
Roseboom et al. 2013, Geach et al. 2016 in preparation). 
In order to enable effective 450 pLm observations, only the 
very best/dryest conditions were used (i.e. r 225 GHz < 0.05), 
and to maximise depth the imaging was undertaken with a 
“daisy” mapping pattern (Bintley et al. 2014). 

The details of the reduction process are described in 
Roseboom et al. (2013), and so only a brief description is 
given here. 

The data were reduced with the SMURF packag^Vl.4.0 
(Chapin et al. 2013) with flux calibration factors (FCFs) of 
606 Jy pW“^ Beam~^ for 450/rm and 556 JypW~^ Beam”^ 
for 850/rm (Dempsey et al. 2013). 

The noise-only maps were constructed by inverting an 
odd half of the ~ 30 min scans and stacking them all 
together. In the science maps the large-scale background 
was removed by applying a high-pass filter above 1.3 Hz 
to the data (equivalent to 120 arcsec given the SCUBA-2 
scan rate). Then a “whitening filter” was applied to sup¬ 
press the noise in the map whereby the Fourier Transform of 
the map is divided by the noise-only map power spectrum, 
normalised by the white-noise level and transformed back 
into real space. The effective point-source response function 
(PRF) was constructed from a Gaussian with a full-width- 
half-maximum (FWHM) of 14.6 arcsec following the same 
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procedure. Finally, the real sources with a signal-to-noise ra¬ 
tio (SNR) of better than five were extracted by convolving 
the whitened map with the above PRF (see §4.2 of Chapin 
et al. 2013). 

The 850 pia image and the sources extracted from it are 
shown in Fig.[^ while the positions and sub-mm photometry 
for the sources are listed in Appendix 0 Table |A1| 

A total of 106 850 pLva sources were found within the 
map with a SNR > 5. The photometry at 450/rm was per¬ 
formed in the same manner,but assuming the PRF at 450 /rm 
to be a Gaussian of FWHM = 8 arcsec. The 450/rm coun¬ 
terparts to the 850 pLva sources were adopted if a 450 ^m- 
selected source was found within 6 arcsec of the 850 pLxa cen¬ 
troid. As seen in Fig. 53 850 /rm sources have 450 /rm 
counterparts with the mean separation of 2.7 ± 0.2 arcsec. 
Otherwise, for the purpose of SED fitting, the 450 fim flux 
density was measured at the 850 p.m position (flags 1 and 0 
in Table Al respectively). 

The completeness of the 850 ^m catalogue was as¬ 
sessed by injecting sources of known flux density into the 
noise-only maps. Overall 10^ objects were used, split into 
10 logarithmically-spaced flux-density bins between 1 and 
60mJy. In total 2000 simulated maps were created and the 
source extraction was performed in the same way as with the 
real maps. The completeness was then assessed by dividing 
the number of extracted sources by the number of sources 
inserted into the noise-only maps, and the results are shown 
in Fig.[g 


2.2 Supporting multi-frequency data 

This first deep S2CLS pointing within the COS¬ 
MOS/UltraVISTA field was chosen to maximise the power 
of the available ancillary multi-wavelength data, in partic¬ 
ular the HST Cosmic Assembly Near-infrared Deep Ex- 
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Figure 3. Completeness of the 850 fim source sample as a func¬ 
tion of flux density based on 2000 simulated maps with source 
extraction performed on 10^ artificially-created objects. 


tragalactic Legacy Survey (CANDELSjj^ imaging (Gro- 
gin et al. 2011). In addition, the optical Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS; Gwyn 2012), 
the Subaru/Suprime-Gam T-band (Taniguchi et al. 2007; 
Furusawa et al., in preparation) and Ultra VISTA near- 
infrared data (McGracken et al. 2012) were used. The cat¬ 
alogues were made by smoothing all the ground-based and 
HST data to the seeing of the UltraVISTA V-band image 
with the Gaussian of FWHM = 0.82 arcsec (for details, see 
Bowler et al. 2012, 2014, 2015). The catalogue was selected 
in the smoothed CANDELS //-band image and photome¬ 
try was measured in 3 arcsec apertures using the dual-mode 
function in SExtractor (Bertin & Arnouts 1996) on all 
other PSF homogenised images. 

The Spitzer IRAC flux densities at 3.6 /rm and 4.5 pm 
were measured from the S-GOSMOS survey (Sanders et al. 
2007), after image deconfusion based on the UltraVISTA Ks- 
band image; using GALFIT (Peng et al. 2002) the iCg-band 
images were modelled, and the corresponding structural pa¬ 
rameters were then applied to both the 3.6 pm and 4.5 pm 
data and the flux-densities allowed to vary until the opti¬ 
mum fit to the IRAG image of each object was achieved 
(after convolution with the appropriate PSFs). The infinite- 
resolution scaled model IRAG images created in this way 
were than smoothed again to match the seeing of the Ul¬ 
traVISTA V-band image, after which the IRAC flrrx densi¬ 
ties were measured within 3 arcsec apertures. For the small 
number of objects selected from the SCUBA-2 map which 
lay outside the area with CANDELS HST imaging (see Fig. 

the As-band UltraVISTA image was used as the primary 
image for near-infrared candidate counterpart selection. 

The 24 pm catalogue was constructed using the MIPS 
24 pm imaging from the S-COSMOS survey (LeFloch et al. 
2009). The source extraction was performed on the publicly- 
available imaging using the STARFINDER idl package (Dio- 
laiti et al. 2000). The resulting catalogue covers ~ 2.1 deg^ 
and reaches the depth of ct ~ 13 ply (for details, see Rose- 
boom et al. 2013). 

For the extraction of far-infrared flrrx densities and lim¬ 


its we used the Herschel (Pilbratt et al. 2010) Multi-tiered 
Extragalactic Survey (HerMES; Oliver et al. 2012) and 
the Photodetector Array Camera and Spectrometer (PACS; 
Poglitsch et al. 2010) Evolutionary Probe (PEP; Lutz et 
al. 2011) data obtained with the Spectral and Photometric 
Imaging Receiver (SPIRE; Griffin et al. 2010) and PACS in¬ 
struments, covering the entire COSMOS field. We utilised 
Herschel maps at 100, 160, 250, 350 and 500 pm with beam 
sizes of 7.39, 11.29, 18.2, 24.9, and 36.3 arcsec, and 5a sen¬ 
sitivities of 7.7, 14.7, 24.0, 27.5, and 30.5 mjy, respectively. 
We obtained the fluxes of each SCUBA-2 source in the fol¬ 
lowing way. We extracted 120-arcsec wide stamps from each 
Herschel map around each SCUBA-2 source and used the 
PACS (100, 160/im) maps to simultaneously fit Gaussians 
with the FWHM of the respective map, centred at all ra¬ 
dio and 24 pm sources located within these cut-outs, and at 
the positions of the SCUBA-2 optical identifications (IDs, 
or just sub-mm positions if no IDs were selected). Then, to 
deconfuse the SPIRE (250, 350 and 500/rm) maps in a simi¬ 
lar way, we used the positions of the 24 pm sources detected 
with PACS (at > 3a), the positions of all radio sources, and 
the SCUBA-2 ID positions. 

Finally the Very Large Array (VLA) COSMOS Deep 
catalogue was used where the additional VLA A-array ob¬ 
servations at 1.4 GHz were obtained and combined with the 
existing data from the VLA-COSMOS Large project (for 
details, see Schinnerer et al. 2010). This catalogue covers ~ 
250 arcmin^ and reaches a sensitivity of ct = 12 /rjy beam~^. 


3 SCUBA-2 SOURCE IDENTIFICATIONS 


In order to find the optical counterparts for sub-mm sources, 
for which positions are measured with relatively large 
beams, a simple closest-match approach is not sufficiently 
accurate. We therefore use the method outlined in Dunlop et 
al. (1989) and Ivison et al. (2007) where we adopt the 2.5 ct 
search radius around the SCUBA-2 position based on the 
signal-to-noise ratio (SNR): = 2.5 x 0.6 x FWHM/SNR, 
where FWHM ~ 15 arcsec. In order to account for sys¬ 
tematic astrometry shifts (caused by pointing inaccuracies 
and/or source blending; e.g. Dunlop et al. 2010) we enforce 
a minimum search radius of 4.5 arcsec. Within this radius we 
calculate the corrected Poisson probability, p, that a given 
counterpart could have been selected by chance. 

For reasons explained below, the VLA 1.4 GHz and 
Spitzer MIPS 21 pm and IRAC 3 pm (with addition of 
3.Q pm) bands were chosen for searching for galaxy coun¬ 
terparts. In the case of the MIPS 24 pm band, the minimum 
search radius was increased to 5 arcsec to account for the 
significant MIPS beam size (~ 6arcsec). The optical/near- 
infrared catalogues were then matched with these coordi¬ 
nates using a search radius of r = 1.5 arcsec and the closest 
match taken to be the optical counterpart. In addition, we 
utilised the Herschel, SCUBA-2 and VLA photometry to 
help isolate likely incorrect identifications (Section |4.2[ ). 

The results of the identification process are summarized 


in Table A2 where the most reliable IDs (p < 0.05) are 
marked in bold, the tentative IDs (0.05 < p < 0.1) are 
marked in italics, and incorrectly identified sources (as dis¬ 


cussed in Section 4.2 1 are marked with asterisks. 


Given the depth of the 850 pm imaging utilised here. 
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Table 1. The radio/infrared/optical identification statistics for the 850/rm S2CLS COSMOS sample used in this work. The number of 
reliably (with a probability of chance association, p < 0.05), tentatively (with 0.05 < p < 0.1) and all (p < 0.1) identified sources are 
given (with the percentage, out of 106 sources, in brackets). The columns give the ID success rate at a given band followed by the overall 
radio/mid-infrared ID success rate, the raw optical ID success rate, and revised optical ID success rate after checking for consistency 
with the long-wavelength photometric redshifts (see Fig. and Section |4.2|l. 



1.4 GHz 

24 pm 

8 pm 

radio/IR 

overall 

optical 
before corr. 

optical 
after corr. 

reliable {p < 0.05) 
tentative (0.05 < p < 0.1) 
all {p < 0.1) 

15 (14%) 

0 (0%) 

15 (14%) 

62 (58%) 
11 (10%) 
73 (69%) 

37 (35%) 
20 (19%) 
57 (54%) 

67 (63%) 
13 (12%) 
80 (75%) 

67 (63%) 

13 (12%) 

80 (75%) 

54 (51%) 

8 (8%) 

62 (58%) 



radial offset / sigma 

Figure 4. A demonstration that the distribution of positional 
offsets between the 850 pm sources and identified galaxy counter¬ 
parts is consistent with statistical expectations. The histograms 
show the cumulative distribution of positional offset divided by 
positional uncertainty, where the positional uncertainty for each 
SCUBA-2-counterpart association is derived by calculating the 
uncertainty in the position of both the 850 pm source and its coun¬ 
terpart (based on the standard formula cr = 0.6 X FWHM/SNR; 
see Section]^, and adding these in quadrature, with an additional 
2arcsec added in quadrature to account for JCMT pointing un¬ 
certainties. The histograms show the distributions for the most 
secure p < 0.05, and less secure p < 0.1 identifications, while 
the curve shows the prediction assuming a Gaussian distribution. 
This plot provides reassurance that, at the flux level of the sources 
considered here, source confusion has not significantly distorted 
the source positions, and source blending is not a significant issue. 


it is important to check that source positions have not 
been significantly distorted by source confusion. We have 
therefore checked that the distribution of positional off¬ 
sets between the 850 pm sources and their adopted multi¬ 
frequency frequency counterparts is as expected, assum¬ 
ing the standard formula for positional uncertainty (i.e. 
a = 0.6 X FWHM/SNR). The results, shown in Fig.pro¬ 
vide reassurance that the vast majority of source positions 
have not been significantly distorted by confusion/blending, 
and that our association process is statistically valid. 


3.1 Radio and 24 pm counterparts 

The 850 pm band is sensitive to the cool dust re-radiating 
energy absorbed from hot, young stars. The radio band also 


traces recent star formation via synchrotron radiation from 
relativistic electrons produced within supernovae (SNe; Con¬ 
don 1992). The 24pm waveband is in turn sensitive to the 
emission from warm dust, and since sub-mm selected galax¬ 
ies are dusty star-forming galaxies, they are also expected 
to be reasonably luminous in this band. There is thus a 
good physical motivation for searching for the counterparts 
of SCUBA-2 sources in the VLA and MIPS imaging. In addi¬ 
tion, the surface density of sources in these wavebands is low 
enough for chance positional coincidences to be rare (given 
a sufficiently small search radius). 

As seen in Table[^ (before the corrections of Section [4.2[ ) 
at 1.4 GHz the ID success rate is only 14% (15 out of 106 
sources, all with p < 0.05) but at 24 pm the success rate is 
69% (73 out of 106, 62 of which have p < 0.05). Combining 
both methods, the successful identification rate is 70% (74 
out of 106, 63 of which have p < 0.05). The striking differ¬ 
ence in these statistics is due to the fact that the S-COSMOS 
24 pm imaging utilised here is relatively deeper than the ra¬ 
dio data currently available in the COSMOS field. 


3.2 8 pm counterparts 

In order to maximise the identification success rate we also 
searched for counterparts in the the S-COSMOS IRAC 8 pm 
imaging. At the redshifts of interest, this waveband traces 
the rest-frame near-infrared light coming from the older, 
mass-dominant stellar populations in galaxies. Given the 
growing evidence that sub-mm galaxies are masssive, it is 
expected they will be more luminous than average in this 
waveband (e.g. Dye et al. 2008; Michalowski et al. 2010; 
Biggs et al. 2011; Wardlow et al. 2011). We found that 57 
of the 106 SCUBA-2 sources (54%) had 8 pm counterparts, 
37 of which have p < 0.05. However, unsurprisingly, several 
of these identifications simply confirmed the identifications 
already secured via the radio and/or 24 pm cross matching, 
and the search for 8 pm counterparts only added 5 new iden¬ 
tifications (2 of which have p < 0.05) to the results described 
in the previous sub-section. 


3.3 Optical counterparts 

In total, therefore, we identified radio/mid-infrared coun¬ 
terparts for 80 of the 106 SCUBA-2 850 pm sources (67 of 
which have p < 0.05; see Table A2|, and hence achieved 
an identification success rate of 75%. The identification suc¬ 
cess rate achieved in each individual waveband is given in 
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Table In addition, we present postage-stamp images for 
all the sources in the online version (Fig. Bl), with all the 
identihcations marked with the appropriate symbols. 

To complete the connection between the SCUBA-2 
sources and their host galaxies, within the area covered by 
the CANDELS HST WFC3/1R imaging (Fig.[3 we matched 
the statistically-significant mid-infrared and radio counter¬ 
parts to the galaxies in the CANDELS iLieo-band imaging 
using a maximum matching radius of 1.5 arcsec. This yielded 
accurate positions for the optical identifications of 60 of 
the SCUBA-2 sources (Table A3l. For those few SCUBA- 
2 sources which lie outside the CANDELS HST imaging, we 
matched the statistically-significant mid-infrared and radio 
counterparts to the galaxies in the TFs-band UltraVISTA 
imaging (using the same maximum matching radius). This 
yielded accurate positions for the optical identifications of 
the remaining 20 sources (Table [A4| ). We note that galax¬ 
ies SC850-37, 46 and 61, even though successfully identified 
in the optical/near-infrared, turned out to be too close to 
a foreground star for reliable photometry (Fig. Bl in the 
online version) and therefore no optical redshifts or stellar 
masses were derived and utilised in the subsequent analysis. 


4 REDSHIFTS 

4.1 Photometric redshifts 

For all the identified sources, the multi-band photometry 
given in Tables |A3| and |A4| was used to derive optical- 
infrared photometric redshifts using a y;^"niinimization 
method (Cirasuolo et al. 2007, 2010) with a template-htting 
code based on the HYPERZ package (Bolzonella et al. 2000). 
To create template galaxy SEDs, the stellar population syn¬ 
thesis models of Bruzual & Chariot (2003) were applied, us¬ 
ing the Chabrier (2003) stellar initial mass function (IMF) 
with a lower and upper mass cut-off of 0.1 and 100 Mq re¬ 
spectively. A range of single-component star-formation his¬ 
tories were explored, as well as double-burst models. Metal- 
licity was fixed at solar, but dust reddening was allowed 
to vary over the range 0 < Ay < 6, assuming the law of 
Calzetti (2000). The HI absorption along the line-of-sight 
was applied according to Madau (1995). The optical-infrared 
photometric redshifts for the 77 optically-identihed sources 
for which photometry could be reliably extracted (i.e. the 80 
identihed sources excluding SC850-37, 46 and 61) are given 
in Table Also given in this table are the optical spectro¬ 
scopic redshifts where available. We note that, in general, 
Zspec and Zp are in excellent agreement, except for the two 
SCUBA-2 sources which are associated with active galac¬ 
tic nuclei (AGN; sources 65 and 72), presumably because 
no AGN template was included in the photometric redshift 
htting procedure. 

In addition, for every SCUBA-2 source we used the 450 
and 850 /rm photometry as well as the Herschel 100, 160, 250, 
350, 500/rm and VLA 1.4 GHz flux densities (or limits) to 
obtain ‘long-wavelength’ photometric redshifts (2lw)- This 
was achieved by fitting the average SED template of sub- 
mm galaxies from Michalowski et al. (2010) to the measured 
flux densities and errors in all 8 of these long-wavelength 
bands (including flux-density measurements corresponding 
to non detections). The resulting ‘long-wavelength’ redshift 
estimates for all 106 sources are also given in Table [AS] 


4.2 Redshift/identification refinement 

Given the statistical nature of the identification process de¬ 
scribed above, there is always a possibility that some identih¬ 
cations are incorrect (as revealed by interferometric follow¬ 
up - e.g. Hodge et al. 2013; Koprowski et al. 2014), and 
indeed, even when the probability of chance coincidence is 
extremely small, it can transpire that the optical counter¬ 
part is not, in fact, the correct galaxy identihcation, but 
is actually an intervening galaxy, gravitationally leasing a 
more distant sub-mm source (e.g. Dunlop et al. 2004). In ei¬ 
ther case, a mis-identihcation will lead to an under-estimate 
of the true redshift of the sub-mm source, and indeed dra¬ 
matic discrepancies between Zp and Zlw can potentially be 
used to isolate mis-identihed sources. 

In Fig. we have therefore plotted zlw versus Zp in 
an attempt to test the consistency of these two indepen¬ 
dent redshift estimators. From this plot it can be seen 
that, for the majority of sources, the two redshift estimates 
are indeed consistent, with the normalized offset in zlw 
(r = [zvff — Zp)/(1 -I- Zp)) displaying a Gaussian distribu¬ 
tion with (7 = 0.14. However, there is an extended pos¬ 
itive tail to this distribution, indicative of the fact that 
a significant subset of the identifications have a value of 
Zp which is much smaller than the (identification indepen¬ 
dent) ‘long-wavelength’ photometric redshift of the SCUBA- 
2 source, zlw- Given the aforementioned potential for mis- 
identihcation (and concomitant redshift under-estimation) 
we have chosen to reject the optical identihcations (and 
hence also Zp) for the sources that lie more than 3cr above the 
1:1 redshift relation (see Fig.j^and caption for details). This 
may lead to the rejection of a few correct identihcations, but 
this is less important than the key aim of removing any sig- 
nihcant redshift biases due to mis-identihcations, and also 
the value of retaining only the most reliable set of identihed 
sources for further study. 

The ehect of this cut is the rejection of 18 of the 80 
optical identihcations derived in Section These rejected 
optical IDs are hagged with asterisks in Table |A2| and ze¬ 
ros in Table |A5| We emphasize that the rejection of these 
low-redshift identihcations does not impact signihcantly on 
the investigation of the physical properties of the sub-mm 
sources at 2 : ~ 1 — 4 pursued further below, because if the 
low-redshift IDs were retained they would not feature in the 
relevant redshift bins, while adoption of the long-wavelength 
redshift for these sources means that we do not include these 
sources in the sample of objects with reliable stellar masses. 
We also stress that only a small subset of these objects are 
likely lenses (5 possible examples are highlighted in Fig. B2 
available in the online version), but while a revised search for 
galaxy counterparts for the other sources might yield alter¬ 
native counterparts with Zp consistent with zlw, we prefer 
not to confuse subsequent analysis by the inclusion of what 
would be inevitably less reliable galaxy identihcations. 

As tabulated in Table with this redshift rehnement, 
the effective optical ID success rate for the most reliable 
{p < 0.05) IDs drops from 63% to 51%, while the overall 
{p < 0.1) ID success rate drops from 75% to 58%. However, 
while this reduces the number of reliably identihed SGUBA- 
2 sources to ~ 50% of the sample, this has the advantage 
or removing the most dubious identihcations. Moreover, we 
stress that we retain redshift information for every one of 



M. Koprowski et al. 




Figure 5. The left-hand panel shows the ‘long-wavelength’ photometric redshifts {zly/) derived for the SCUBA-2 sources plotted against 
the optical-infrared photometric redshifts {zp) of the optical identifications (see Section |4.2[ |. The central blue solid line shows the 1:1 
relation. As illustrated in the right-hand panel, the sources lying below the 1:1 relation display a distribution of normalized redshift offsets 
(i.e. r = (zlw ~ 2 :p)/(l -|- Zp)) which is approximately Gaussian with a = 0.14. The positive side of this distribution is also reasonably 
well fitted by this same Gaussian, but there is a long positive tail, indicative of the fact that a significant subset of the identifications 
have a value of Zp which is much smaller than the (identification independent) ‘long-wavelength’ photometric redshift of the SCUBA-2 
source Given the potential for mis-identification (e.g. through galaxy-galaxy gravitional lensing) we view such discrepancies as 

evidence that Zp, or more likely the galaxy identification itself, is in error. The upper and lower blue solid lines in the left-hand panel 
show the di3(j limits of the Gaussian distribution, and so we choose to reject the optical identifications (and hence also Zp) for the sources 
that lie above the Scr limit (red dots). This same 3cr limit is shown by the black vertical line in the right-hand panel. 
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Figure 6. Left-hand panel: The redshift distribution of our full 106-source S2CLS 850//m sample in the COSMOS field (Table [A^ . 
The black area shows the distribution for the robust optical identifications with spectroscopic or optical-infrared photometric redshifts, 
which has a mean redshift of z = 1.97 di 0.09. The blue line depicts the redshift distribution of the enlarged sample which results from 
adding the SCUBA-2 sources which lack robust optical identifications, but which have reliable ‘long-wavelength’ redshifts (z^wr)- Finally, 
the additional histrogram maked by the green arrows indicates the objects for which only lower-limits on redshift could be derived from 
the long-wavelength photometry. Adopting these lower limits, the mean redshift for the whole sample is z = 2.38 di 0.09. Right-hand 
panel: The redshift distribution for the whole S2CLS COSMOS sample overlaid with the redshift distributions derived by Chapman et 
al. (2005) {z = 2.00 di 0.09), and for the robust galaxy identifications in the AzTEC/SHADES survey presented by Michalowski et al. 
(2012a) (z = 2.00 di 0.10). In addition we plot the redshift distribution of the sample of luminous (sub-)mm sources in the COSMOS 
field presented by Koprowski et al. (2014) (z = 3.53 di 0.19). 



the 106 SCUBA-2 sources, in the form of zlw if neither 
Zspec nor a reliable value for Zp are available. 

In Fig. B2 (available in the online version) we present 
12 X 12 arcsec near-infrared postage-stamp images for every 
source, with the positions of all the IDs marked. In this 
figure we give the source name in red if the optical ID was 


in fact subsequently rejected in the light of zlw- It can be 
seen from this figure that at least some of these incorrect 
identifications are indeed due to galaxy-galaxy lensing (see 
figure caption for details). 
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Figure 7. 850 fim flux density as a function of redshift for the whole SCUBA-2 850-/i.m sample used in this study. The open symbols 
show the 80 sources which possess optical IDs and the filled circles show the 26 sources that lack optical IDs, with arrows signifying lower 
limits on ‘long-wavelength’ redshifts. All the objects are colour-coded according to their 850 fim signal-to-noise-ratio. It can be seen that 
the galaxies that lack optical associations span a similar range in SNR to the optically identified sources. This suggests that the failure to 
uncover optical counterparts for these 26 sources is not due to them being spurious sub-mm detections, but instead suggests that either 
they lie at somewhat higher redshifts than the typical identifled source, and/or that these unidentified SCUBA-2 sources are blends of 
individual fainter sources whose counterparts lie below the radio/24 pm detection threshold. 


Table 2. Five sources in our sample that have been the subject of previous detailed study. Four of them (the ones with Zprev) were 
previously followed up with the mm/sub-mm interferometry. The columns show respectively our ID, the source name from previous work 
(the full previous ID name for the AzTEC source is AzTEC J100025.23+022608.0), interferometric RA and Dec (or, in the case of the 
AzTEC source, single-dish coordinates), the separation between the interferometric position (except AzTEC) and the optical ID found in 
this work, our final redshifts (from Table [A^ , the redshift estimates from previous studies, and finally references. Where two references 
are given, the first one refers to the coordinates and the second to the previous redshift estimate (i.e. the source of Zprev)- 


ID 

IDother 

RA 

Dec 

separation 

z 

Zprev 

References 

SC850- 


/deg 

/deg 

/arcsec 




1 

MMl 

150.0650 

2.2636 

2.62 

3.30+?l 

Q 1 n+0-50 

Aravena et al. (2010) 

6 

COSLA-35 

150.0985 

2.3653 

0.13 

2.50+?° 

Q 1 fi+O-24 

'^•49_o.26 

Smolcic et al. (2012); Koprowski et al. (2014) 

14 

COSLA-8 

150.1064 

2.2523 

2.76 

2.18++ 

i.yu_0.22 

Smolcic et al. (2012); Koprowski et al. (2014) 

29 

AzTEC 

150.1051 

2.4356 

2.03 

2 41 +0-24 


Scott et al. (2008) 

31 

COSLA-38 

150.0525 

2.2456 

0.27 

9 /17+0-0® 
-0.12 

2 

^•^^-0.11 

Smolcic et al. (2012) 
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4.3 Redshift distribution 

The differential redshift distribution for our SCUBA-2 
galaxy sample is presented in Fig. In the left-hand panel 
the black area depicts the redshift distribution for the 
sources with reliable optical IDs (and hence Zspec or Zp), 
while the histrogram indicated in blue includes the addi¬ 
tional unidentified SCUBA-2 sources with meaningful mea¬ 
surements of zlw- Finally, the green histogram containing 
the green arrows indicates the impact of also including those 
sources for which only lower limits on their estimated red- 
shifts could be derived from the long-wavelength photom¬ 
etry. The mean and median redshifts for the whole sam¬ 
ple are 2 = 2.38 ± 0.09 (strictly speaking, a lower limit) 
and Zmed = 2.21 ± 0.06 respectively whereas, for the con¬ 
firmed optical IDs with optical spectroscopic/photometric 
redshifts the corresponding numbers are 2 = 1.97 ± 0.09 
and Zmed = 1.96 ± 0.07. This shows that, as expected, the 
radio/infrared identification process biases the mean red¬ 
shift towards lower redshifts, but in this case only by about 
~ 10% in redshift. In addition, to make sure that our uniden¬ 
tified sources are in fact not spurious, which would manifest 
itself as them having low SNR values, we also plot in Fig.[^ 
the 850 fim flux as a function of redshift for the whole sample 
used here, colour-coded according to SNR. It can be clearly 
seen that the unidentified sources exhibit a wide range of 
SNRs and hence are most likely real. 

In the right-hand panel of Fig. we compare the red¬ 
shift distribution of the deep 850 pm selected sample stud¬ 
ied here with example redshift distributions from previous 
studies. Although our sample is somewhat deeper/fainter 
than the sub-mm samples studied previously by Chapman 
et al. (2005) and Michalowski et al. (2012a), the redshift 
distributions displayed by the optically-identified subset of 
sources from each study are remarkably consistent; we find 
2 = 1.97 ± 0.09, while Chapman et al. (2005) reported 
2 = 2.00 ± 0.09, and Michalowski et al. (2012a) reported 
2 = 2.00 ± 0 . 10 . 

While inclusion of our adopted values of zlw for our 
unidentified sources moves the mean redshift up to at least 
2 ~ 2.4, it is clear that the redshift distribution found here 
cannot be consistent with that found by Koprowski et al. 
(2014) for the subset of very bright sub-mm/mm sources in 
the COSMOS field (see also Smolcic et al. 2012), for which 
2 = 3.53±0.19. This is not due to any obvious inconsistency 
in redshift estimation techniques, as can be seen from Table 
[^(discussed further below), and indeed the analysis methods 
used here are near identical to those employed by Koprowski 
et al. (2014). Rather, as discussed in Koprowski et al. (2014), 
there must either be a trend for the most luminous sub- 
mm/mm sources (i.e. Sgso^m > 8 mJy) to lie at significantly 
higher redshifts than the more typical sources studied here, 
or the COSMOS bright source sample of Scott et al. (2008) 
imaged by Younger et al. (2007, 2009) and Smolcic et al. 
( 2012 ) must be unusually dominated by a high-redshift over¬ 
density in the COSMOS field. 


4.4 Previous literature associations 

Five of the sub-mm sources in our SCUBA-2 sample have 
been previously studied in some detail, and so, in Table[^ we 
compare our ID positions and redshifts with the pre-existing 


information. Four of these bright sources were previously the 
subject of interferometric mm/sub-mm observations, yield¬ 
ing robust optical identifications and photometric redshifts 
in good agreement with our results. The source separation 
for SC850-29 (2.03 arcsec) is perfectly plausible since this is 
the separation between the original AzTEC single-dish coor¬ 
dinate and our chosen ID. The small separations between the 
positions of our adopted IDs for SC850-6 and 31 and their 
mm/sub-mm interferometric centroids confirm the reliabil¬ 
ity of our ID selection. For SC850-1 the rather large source 
separation of 2.62 arcsec supports our rejection of the opti¬ 
cal ID for this source. Finally, the rather large separation for 
SC850-14 clearly casts doubt on our adopted ID, but in this 
case Zp is very similar to zlw (which, of course, is why we 
did not reject the ID) and so the final redshift distribution 
is unaffected by whether or not the ID is correct. 


5 PHYSICAL PROPERTIES 


5.1 Stellar masses and star-formation rates 


For the 58 SCUBA-2 sources for which we have secure op¬ 
tical identifications±redshifts (after the sample refinement 
discussed in Section 4.2 1 we were able to use the results 
of the SED fitting (used to determine Zp) to obtain an es¬ 
timate of the stellar mass, M^,, for each galaxy. The de¬ 
rived stellar masses were based on the models of Bruzual 
& Chariot (2003) assuming double-burst star-formation his¬ 
tories (see Michalowski et al. 2012b), and we assumed a 
Chabrier (2003) IMF. 


We were also able to estimate the star-formation rate, 
SFR, for each of these sources by using the average long- 
wavelength SED of the sub-mm galaxies from Michalowski et 
al. (2010), applied to the 850/rm flux-density of each source 
at the relevant photometric redshift, to estimate the far- 
infrared luminosity of each source. 

The resulting SFRs are plotted against M* in Fig.[^ In 
the main plot, for clarity we have confined attention to the 
sources with Zp > 1.5 because, as shown in the inset plot, due 
to the impact of the negative K-correction at 850 fim, at 2 > 

1.5 the flux-density limit of the current sample essentially 
equates to SFR ~ 100 MQyr”'^ at all higher redshifts. In this 
plot we also show the position of the ‘Main-Sequence’ (MS) 
of star-forming galaxies, as deduced at 2 ~ 2.5 by Elbaz et 
al. (2011), and at 2 > 1.5 by Rodighiero et al. (2011). The 
sensitivity of our deep SCUBA-2 sample to values of SFR 
as low as lOOM 0 yr“^ means that, for objects with stellar 
masses M* > 7 x 10^° M©, we are able for the first time to 
properly compare the positions of sub-mm selected galaxies 
on the SFR:M* plane with the MS in an unbiased manner. 


5.2 Specific star-formation rates 

In Fig. 1^ we collapse the information shown in Fig. into 
distributions of specific SFR. The black histogram shows the 
distribution of sSFR for the whole robustly-identified sample 
of SCUBA-2 sources at 2 > 1.5, but this can be subdivided 
by mass into the sub-sample with M* > 7 x 10^° Mq (blue 
histogram) and the complementary sub-sample of sources 
with M* < 7 X 10^° M 0 (red histogram). Refering back to 
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Figure 8. The star-formation rate (SFR) as a function of stellar mass (M*) for the robustly identified SCUBA-2 sources with 2 : > 1.5. 
As can be seen from the inset plot, due to the impact of the negative K-correction at 850/im, for 2 > 1.5 the flux-density limit of 
the current sample essentially equates to SFR ~ lOOMQyr”^. The black solid line in the main plot shows the position of the so-called 
‘Main-Sequence’ (MS) of star-forming galaxies at 2 ~ 2.5 as deduced by Elbaz et al. (2011), while the black dashed line depicts the MS 
at 2 > 1.5 as given by Rodighiero et al. (2011). The sensitivity of our deep SCUBA-2 sample to values of SFR as low as lOOMQyr”^ 
means that, for objects with stellar masses M* > 7 x lO^^MQyr"^ (i-®-? inside the cyan rectangle), we are able for the first time to 
properly compare the positions of sub-mm selected galaxies on the SFR:M* plane with the MS in an unbiased manner. As shown in 
Fig. H we find that, confining our attention to M* > 7 x 10^*^ M 0 yr the SCUBA-2 sources display a Gaussian distribution in specific 
SFR peaking at sSFR = 2.25 di 0.19 Gyr“^ (corresponding to the main sequence shown here by the blue solid line), demonstrating that 
the SCUBA-2 sources lie on the high-mass end of the normal star-forming MS at 2 ~ 2. 


Fig.§ it can be seen that, at lower stellar masses, the mea¬ 
surement of sSFR is inevitably biased high by the effective 
SFR limit > lOOMQyr"^, and so it is difficult to tell if 
these SCUBA-2 sources genuinely lie above the MS, or if 
we are simply sampling the high-sSFR tail of the distribu¬ 
tion around the MS. However, at M* > 7 x 10^® Mq it is 
clear that the SFR limit would not produce a significantly 
biased sampling of the distribution of galaxies on the MS. In 
essence, because of the depth of the SCUBA-2 imaging, for 
sub-mm selected galaxies with M* > 7 x 10^'^ M© we should 
now be able to perform the first unbiased estimate of their 
sSFR at 2 ~ 1.5 — 3. 

In fact, for the high-mass sub-sample, in which SFR 
is not biased by the effective flux-density limit of the 
deep SCUBA-2 survey, the distribution of sSFR resem¬ 
bles closely a Gaussian peaked at sSFR = 2.25 Gyr“^ with 
a = 0.89 Gyr“^. This Gaussian fit is shown by the green 
curve in Fig. and is completely consistent with the nor¬ 
malization and scatter (~ 0.25 dex) in the MS reported by 
Rodighiero et al. (2011). 

Finally, to check whether we could be biased towards 
high-mass (and hence low sSFR) objects at high redshift, as 
a consequence of the flux-density limits of our optical/near- 
infrared catalogues, we plot the near-infrared (CANDELS 


H-band and UltraVISTA ifs-band) absolute magnitudes of 
our source IDs against redshift in Fig. |10| The measured 
values are generally not close to the detection limits of our 
catalogues and therefore we conclude that the sample is not 
biased against high sSFRs at high redshifts on account of 
an inability to detect low-mass galaxies. 

We conclude, therefore, within the stellar mass range 
where we are able to sample the distribution of sSFR in an 
unbiased way, the sub-mm sources uncovered from this deep 
SCUBA-2 850 fim image, display exactly the mean sSFR and 
scatter expected from galaxies lying on the high-mass end 
of the star-forming main-sequence at z ~ 2. 


5.3 The ‘main sequence’ and its evolution 

Given that the SCUBA-2 sources seem to, in effect, define 
the high-mass end of the star-forming main sequence (MS) 
of galaxies over the redshift range probed by our sample 
(i.e. 1.5 < z < 3) it is of interest to explore how the inferred 
normalization and slope of the MS as derived here compares 
to that derived from other independent studies based on very 
different selection techniques over a wide range of redshifts. 

Thus, in Fig. we divide our (high-mass) sample into 
three redshift bins to place the inferred evolution of sSFR 
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Figure 9. The distribution of specific star-formation rate, sSFR, 
as derived from the values of SFR and M* plotted in Fig. The 
black histogram shows the distribution for the whole robustly- 
identified sample of SCUBA-2 sources at 2 > 1.5, but this can be 
subdivided by mass into the sub-sample with M* > 7 x 10^® Mq 
(blue histogram) and the complementary sub-sample of sources 
with M* < 7 X 10^^ Mq (red histogram). It can be seen that, for 
the high-mass sample, in which SFR is not biased by the effective 
flux-density limit of the deep SCUBA-2 survey, the distribution 
resembles closely a Gaussian peaked at sSFR = 2.25 Gyr“^ with 
cr = 0.89 Gyr“^, as shown by the green curve. This demonstrates 
that, where their distribution on the SFR:M* plane can now fi¬ 
nally be probed in an unbiased manner, the SCUBA-2 galaxies 
lie on the MS of star-forming galaxies at 2 ~ 2. 



within the wider context of studies spanning virtually all of 
cosmic time (i.e. 0 < 2 < 8 ). 

The first obvious striking feature of Fig. EH is that our 
new determination of average sSFR over the redshift range 
1.5 < a < 3 follows very closely the trend defined by the orig¬ 
inal studies of the MS undertaken by Noeske et al. (2007) 
and Daddi et al. (2007). Since such studies were based on 
very different samples, sampling lower stellar masses, this 
result also implies that we Hnd no evidence for a high-mass 
turnover in the MS at these redshifts (i.e. a decline in sSFR, 
or change in the slope of the MS above some characteristic 
mass). Evidence for a decline in the slope of the MS above 
a stellar mass log(M*/M 0 ) ~ 10.5 has been presented by 
several authors (e.g. Whitaker et al. 2014; Tasca et al. 2015) 
but these results are based on optical/near-infrared studies, 
and suffer from two problems. First, as recently discussed by 
Johnston et al. (2015), the results of optically-based studies 
depend crucially on how one selects star-forming galaxies, 
and colour selection can yield an apparent turn-over in the 
MS at high masses simply due to increased contamination 
from passive galaxies/bulges (see also Renzini & Peng 2015; 
Whitaker et al. 2015). Second, and more important, at the 
high SFRs of interest here, it is well known that SED fit¬ 
ting to optical-infrared data struggles to capture the total 
star-formation rate because the vast majority of the star- 
formation activity in high-mass galaxies is deeply obscured. 
It is therefore interesting that other recent studies of the MS 
based on far-infrared/sub-mm data also find no evidence for 
a high mass turnover in the MS at high-redshift; for example 
Schreiber et al. (2015), from their Herschel stacking study 
of the MS, report that any evidence for a flattening of the 


Figure 10. Absolute magnitude versus redshift for the secure 
SCUBA-2 near-infrared IDs as derived for the sources lying within 
the HST imaging (i.e. CANDELS i7-band neasurements; upper 
panel) and those lying outside the HST imaging (i.e. UltraVISTA 
iiT^-band measurements; lower panel). The blue and green solid 
curves represent the detection limits of our CANDELS 77-band 
and the UltraVISTA Tfs-band selected catalogues respectively. 
It can be seen that virtually all the sources lie well above the 
detection limits, indicating that this our sample is not vulnerable 
to serious biases or incompleteness in stellar mass. 

MS above log(M*/M 0 ) ~ 10.5 becomes less prominent with 
increasing redshift and vanishes by z ~ 2 . 

As is clear from Fig.[^ the present study does not pro¬ 
vide sufficient dynamic range to enable a new measurement 
of the precise value and redshift evolution of the slope of 
the MS (see Speagle et al. 2014 for results from a compila¬ 
tion of 25 studies). Nevertheless, the advantages of sub-mm 
selection for an unbiased study of the high-mass end of the 
MS are clear (i.e. no contamination from passive galaxies, 
and a complete census of dust-enshrouded star formation), 
and our results show that the slope of the MS must remain 
close to unity up to stellar masses M* ~ 2 x 10 ^^ M 0 at 
2 ; ~ 2 — 3. We note that it is sometimes claimed that stud¬ 
ies of the MS based on far-IR or sub-mm selected samples 
yield vastly different determinations of the SFR-M* rela¬ 
tion from the MS (e.g. Rodighiero et al. 2014), but it needs 
to be understood that this is because previous studies based 
on such samples did not reach sufficient sensitivity in SFR 
(for individual objects) to properly sample the MS at high 
redshift. As emphasized in Section [5.2[ and in Fig. even 
the deepest ever 850/rm survey analysed here only enables 
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Figure 11. Mean sSFR as a function of the redshift. The results of this work (inset plot with the mean values represented by black 
points with thick error bars) calculated using the unbiased sample (from inside the cyan rectangle in Fig.j^ are shown by the magenta 
diamonds. It can be seen that the evolution of characteristic sSFR (= to the normalization of the MS) inferred from the SCUBA-2 
galaxies is in excellent accord with the results from many other studies. 


US to properly explore the MS at the very highest masses, 
due to the effective SFR sensitivity limit; clearly the sources 
detected in the present study at lower masses are outliers 
from the MS, and can only provide indirect information of 
the scatter in the MS at masses of a few xIO^^Mq, rather 
than its normalization. 

Finally, looking to higher redshifts, Fig. |11| shows that 
the present study does not provide useful information on 
characteristic sSFR beyond 2 ~ 3, but also demonstrates 
that the trend indicated here over 1.5 < 2 < 3 extends nat¬ 
urally out to our previous determination of sSFR in very 
high-redshift sub-mm/mm galaxies at 2 ~ 5 (Koprowski 
et al. 2014). There is currently considerable debate over 
the normalization of the MS at 2 ~ 4, due in large part 
to uncertainty over the impact of nebular emission lines 
on the estimation of stellar masses (see e.g. Stark et al. 
2013; Smit et al. 2014). However, the sub-mm studies of 
high-mass star-forming galaxies are clearly consistent with 
the results of several existing studies (e.g. Steinhardt et al. 
2014), and (despite their supposedly extreme star-formation 
rates) sub-mm-selected galaxies provide additional support 
for the presence of a ‘knee’ in the evolution of sSFR around 
2 ~ 2 (as originally suggested by the results of Gonzalez et 
al. 2010, 2012). The ability of theoretical models of galaxy 
formation to reproduce this transition remains the subject 
of continued debate, with smooth cold accretion onto dark 
matter halos leading to expectations that sSFR should rise 
oc (1 + 2 )^'® (Dekel et al. 2009, 2013; Facher-Giguere et 


al. 2011; Rodn'guez-Puebla et al. 2016), and a range of hy- 
drodynamical and semi-analytic models of galaxy formation 
yielding predictions of characteristic sSFR at 2 ~ 2 that fall 
short of the results shown in Fig. 11 by a factor of 2 — 6 (see 


discussion in Johnston et al. 2015, and references therein). 
However, in Fig. we show that when the redshift axis 
is re-cast in terms of cosmic time, there is really no obvi¬ 
ous feature in the evolution of characteristic sSFR. Rather, 
the challenge for theoretical models is to reproduce the ap¬ 
parently simple fact that logj^psSFR is a linear function of 
the age of the Universe, at least out to the highest redshifts 
probed to date. 


6 SUMMARY 

We have investigated the multi-wavelength properties of the 
galaxies selected from the deepest 850-/rm survey under¬ 
taken to date with SGUBA-2 on the JCMT. This deep 850- 
/rm imaging was taken in parallel with deep 450 /rm imaging 
in the very best observing conditions as part of the SGUBA- 
2 Gosmology Legacy Survey. A total of 106 sources (>5cr) 
were uncovered at 850/im from an area of ~ 150arcmin^ in 
the centre of the GOSMOS/UltraVfSTA/CANDELS field, 
imaged to a typical depth of ugso — 0.25 mjy. Aided by 
radio, mid-IR, and 450-/rm positional information, we es¬ 
tablished statistically-robust galaxy counterparts for 80 of 
these sources (~ 75%). 
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Figure 12. Mean sSFR as a function of cosmic time. Data/symbols are as in Fig. The straight-line fit to the data has the form 
logio (sSFR/Gyr-1) = -0.12(t/Gyr) + 0.71. 


By combining the optical-infrared photometric 
redshifts, Zp, of these galaxies with independent ‘long- 
wavelength’ estimates of redshift, zlw (based on 
i/ersc/iel/SCUBA-2/VLA photometry), we have been 
able to refine the list of robust galaxy identifications. This 
approach has also enabled us to complete the redshift 
content of the whole sample, yielding z = 2.38 ± 0.09, a 
mean redshift comparable with that derived from all but 
the brightest previous sub-mm samples. 

Because our new deep 850-/im selected galaxy sam¬ 
ple reaches flux densities equivalent to star-formation rates 
SFR ~ 100 M 0 yr”^, we have been able to confirm that sub- 
mm galaxies form the high-mass end of the ‘main sequence’ 
(MS) of star-forming galaxies at z > 1.5 (with a mean spe¬ 
cific SFR of sSFR = 2.25 ± 0.19Gyr“^ at 2 ~ 2.5). Our 
results are consistent with no significant flattening of the 
MS towards high stellar masses at these redshifts (i.e. SFR 
continues oc M*), suggesting that reports of such flatten¬ 
ing are based on contamination by passive galaxies/bulges, 
and/or under-estimates of dust-enshrouded star-formation 
activity in massive star-forming galaxies. However, our find¬ 
ings contribute to the growing evidence that average sSFR 
rises only slowly at high redshift, from sSFR ~ 2Gyr“^ at 
2 ; ~ 2 to sSFR ~ 4 Gyr“^ at 2 ~ 5. These results are consis¬ 
tent with a rather simple evolution of global characteristic 
sSFR, in which log^gsSFR is a linear function of the age of 
the Universe, at least out to the highest redshifts probed to 
date. 
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APPENDIX A: DATA TABLES 

In this appendix we provide tables detailing: i) the sub-mm 
properties of the deep 106-source 850/rm-selected SCUBA-2 
sample utilised in this study, ii) the results of the galaxy 
counterpart identification process, iii) the optical-infrared 
photometry for the galaxy identifications, and iv) the esti¬ 
mated redshifts and derived physical properties of the sub- 
mm galaxies. 
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Table Al. The basic properties of the 106-source SCUBA-2 850 ^m-selected sample. The penultimate column gives the SCUBA-2 colour 
where, if the significance of the 450 //m detection is less than 2(j, the SCUBA-2 colour is based on a 450 /i.m limit with ^450 < ‘S '450 + 2(j. 
The flag given in the flnal column indicates whether the 450//m flux density was taken from 450 ^m catalogue (1) or simply measured 
at the 850 ^m position (0); the latter measurement was adopted if no 450/^m-selected source with ^450 > 4(j was found within 6arcsec 
of the 850/im source position). 


ID 

RAgso 

/deg 

DECgso 

/deg 

Ssso 

/mjy 

ASgso 

/mJy 

SNRgso 

'S '450 

/mJy 

AS 450 

/mJy 

SNR 450 

'S850/‘S450 

flag 

1 

150.06518 

2.26412 

15.64 

0.38 

41.69 

26.99 

2.38 

11.34 

0.58 

1 

2 

150.09985 

2.29772 

10.20 

0.28 

36.82 

17.74 

1.77 

10.00 

0.58 

1 

3 

150.10079 

2.33499 

7.33 

0.23 

32.02 

10.32 

1.41 

7.34 

0.71 

1 

4 

150.10549 

2.31327 

7.79 

0.24 

31.96 

23.42 

1.53 

15.35 

0.33 

1 

5 

150.14320 

2.35607 

7.88 

0.26 

29.98 

19.71 

1.54 

12.78 

0.40 

1 

6 

150.09833 

2.36568 

8.20 

0.28 

29.20 

22.81 

1.80 

12.71 

0.36 

1 

7 

150.09847 

2.32162 

7.04 

0.25 

28.44 

16.66 

1.53 

10.88 

0.42 

1 

8 

150.09820 

2.26061 

6.44 

0.33 

19.32 

14.89 

2.13 

6.98 

0.43 

1 

9 

150.07809 

2.28168 

5.88 

0.32 

18.56 

15.45 

2.03 

7.62 

0.38 

1 

10 

150.15390 

2.32833 

4.75 

0.26 

18.17 

11.15 

1.55 

7.17 

0.43 

1 

11 

150.04264 

2.37371 

7.34 

0.41 

17.85 

23.66 

2.87 

8.23 

0.31 

1 

12 

150.10996 

2.25832 

5.54 

0.34 

16.55 

8.91 

2.13 

4.18 

0.62 

1 

13 

150.08512 

2.29050 

4.87 

0.30 

16.25 

12.79 

1.97 

6.50 

0.38 

1 

14 

150.10692 

2.25218 

5.83 

0.36 

16.03 

23.81 

2.38 

9.99 

0.24 

1 

15 

150.11717 

2.33026 

3.41 

0.21 

15.95 

6.53 

1.31 

4.99 

0.52 

1 

16 

150.05633 

2.37363 

5.42 

0.35 

15.53 

24.31 

2.39 

10.17 

0.22 

1 

17 

150.20799 

2.38297 

7.34 

0.47 

15.50 

15.66 

2.68 

5.83 

0.47 

1 

18 

150.16393 

2.37274 

6.13 

0.40 

15.43 

31.22 

1.96 

15.90 

0.20 

1 

19 

150.11258 

2.37633 

4.35 

0.29 

15.14 

9.91 

1.84 

5.37 

0.44 

1 

20 

150.15024 

2.36457 

4.55 

0.31 

14.46 

13.31 

1.76 

7.56 

0.34 

1 

21 

150.09873 

2.31118 

3.68 

0.26 

13.89 

9.08 

1.62 

5.61 

0.41 

1 

22 

150.05727 

2.29352 

4.60 

0.33 

13.88 

12.83 

2.19 

5.87 

0.36 

1 

23 

150.12283 

2.36081 

3.16 

0.24 

13.36 

10.41 

1.50 

6.94 

0.30 

1 

24 

150.10937 

2.29455 

3.43 

0.27 

12.58 

11.88 

1.74 

6.84 

0.29 

1 

25 

150.03791 

2.34079 

4.56 

0.36 

12.49 

9.20 

2.48 

3.71 

0.50 

0 

26 

150.08011 

2.34091 

3.27 

0.27 

11.97 

10.70 

1.77 

6.06 

0.31 

1 

27 

150.17416 

2.35283 

4.07 

0.34 

11.83 

8.82 

2.01 

4.38 

0.46 

1 

28 

150.12169 

2.34175 

2.48 

0.21 

11.59 

10.56 

1.30 

8.13 

0.23 

1 

29 

150.10535 

2.43531 

6.47 

0.57 

11.31 

18.98 

3.59 

5.28 

0.34 

1 

30 

150.14489 

2.37645 

3.37 

0.32 

10.54 

8.44 

1.78 

4.75 

0.40 

1 

31 

150.05250 

2.24477 

7.85 

0.76 

10.40 

14.21 

5.63 

2.52 

0.55 

0 

32 

150.06641 

2.41264 

4.72 

0.46 

10.29 

6.49 

3.18 

2.04 

0.73 

0 

33 

150.04153 

2.28039 

4.01 

0.42 

9.53 

6.73 

2.60 

2.59 

0.60 

0 

34 

150.13514 

2.39948 

3.03 

0.32 

9.39 

11.12 

1.94 

5.73 

0.27 

1 

35 

150.16742 

2.29950 

3.29 

0.35 

9.36 

11.12 

1.91 

5.84 

0.30 

1 

36 

150.08208 

2.41590 

3.95 

0.43 

9.11 

10.88 

2.91 

3.74 

0.36 

0 

37 

150.06812 

2.27618 

3.06 

0.34 

9.08 

11.70 

2.08 

5.62 

0.26 

1 

38 

150.07620 

2.38036 

3.14 

0.35 

8.87 

12.41 

2.27 

5.46 

0.25 

1 

39 

150.09322 

2.24697 

3.69 

0.43 

8.63 

11.67 

2.83 

4.12 

0.32 

1 

40 

150.10570 

2.32638 

1.94 

0.23 

8.52 

7.87 

1.38 

5.72 

0.25 

1 

41 

150.12888 

2.28474 

2.47 

0.29 

8.48 

0.96 

1.91 

0.50 

> 0.51 

0 

42 

150.02819 

2.34702 

3.80 

0.45 

8.36 

0.76 

2.97 

0.26 

> 0.57 

0 

43 

150.17214 

2.24149 

4.97 

0.59 

8.35 

3.39 

3.87 

0.88 

> 0.45 

0 

44 

150.13663 

2.23305 

4.66 

0.57 

8.23 

2.28 

3.48 

0.66 

> 0.50 

0 

45 

150.12744 

2.38798 

2.52 

0.31 

8.21 

4.50 

1.87 

2.41 

0.56 

0 

46 

150.10606 

2.42844 

3.94 

0.48 

8.14 

14.48 

3.09 

4.69 

0.27 

1 

47 

150.04863 

2.25278 

4.95 

0.62 

8.04 

3.14 

4.64 

0.68 

> 0.40 

0 

48 

150.02227 

2.28899 

5.05 

0.63 

8.00 

16.71 

3.63 

4.60 

0.30 

1 

49 

150.15725 

2.35741 

2.58 

0.33 

7.83 

10.44 

1.84 

5.67 

0.25 

1 

50 

150.08103 

2.36298 

2.39 

0.31 

7.79 

0.75 

1.95 

0.38 

> 0.51 

0 

51 

150.03551 

2.28537 

3.56 

0.46 

7.76 

1.63 

2.76 

0.59 

> 0.50 

0 

52 

150.03693 

2.31959 

2.85 

0.37 

7.72 

1.68 

2.47 

0.68 

> 0.43 

0 

53 

150.18780 

2.32296 

2.54 

0.33 

7.64 

15.59 

1.84 

8.46 

0.16 

1 
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Table Al. (continued). 


ID 

RAgso 

/deg 

DECgso 

/deg 

S'sso 

/mjy 

ASgso 

/mJy 

SNRgso 

>5450 

/mJy 

AS 450 

/mJy 

SNR 450 

>S’850/>S450 

flag 

54 

150.04307 

2.29982 

2.87 

0.36 

7.63 

10.50 

2.36 

4.45 

0.27 

1 

55 

150.13442 

2.37059 

2.06 

0.27 

7.55 

10.13 

1.68 

6.04 

0.20 

1 

56 

150.05005 

2.38574 

3.09 

0.41 

7.50 

7.69 

2.95 

2.61 

0.40 

0 

57 

150.15614 

2.41984 

3.38 

0.45 

7.30 

2.77 

2.72 

1.02 

> 0.41 

0 

58 

150.10709 

2.34444 

1.62 

0.22 

7.24 

4.98 

1.36 

3.66 

0.32 

0 

59 

150.18368 

2.38879 

2.83 

0.39 

7.22 

10.40 

2.14 

4.86 

0.27 

1 

60 

150.19199 

2.27300 

3.25 

0.46 

7.12 

0.68 

2.95 

0.23 

> 0.49 

0 

61 

150.05419 

2.39615 

3.06 

0.43 

7.10 

8.03 

3.08 

2.61 

0.38 

0 

62 

150.16689 

2.23608 

4.61 

0.65 

7.08 

8.03 

4.13 

1.95 

> 0.28 

0 

63 

150.07608 

2.39821 

2.70 

0.39 

6.85 

0.89 

2.62 

0.34 

> 0.44 

0 

64 

150.13004 

2.31505 

1.59 

0.23 

6.82 

8.63 

1.45 

5.95 

0.18 

1 

65 

150.09167 

2.39837 

2.71 

0.40 

6.78 

11.02 

2.49 

4.42 

0.25 

1 

66 

150.17480 

2.40168 

2.70 

0.40 

6.76 

9.29 

2.18 

4.27 

0.29 

1 

67 

150.11157 

2.40409 

2.38 

0.35 

6.73 

3.22 

2.27 

1.42 

> 0.31 

0 

68 

150.13019 

2.25338 

2.37 

0.36 

6.68 

7.01 

2.29 

3.06 

0.34 

0 

69 

150.15507 

2.24389 

3.10 

0.47 

6.63 

-2.98 

3.02 

-0.99 

> 0.51 

0 

70 

150.02490 

2.29668 

3.43 

0.52 

6.63 

3.49 

3.16 

1.11 

> 0.35 

0 

71 

150.07211 

2.23837 

4.44 

0.67 

6.62 

-6.58 

4.78 

-1.38 

> 0.46 

0 

72 

150.06512 

2.32922 

1.93 

0.29 

6.60 

8.04 

2.09 

3.84 

0.24 

0 

73 

150.20910 

2.35567 

2.82 

0.43 

6.60 

18.21 

2.40 

7.57 

0.15 

1 

74 

150.07115 

2.30605 

2.07 

0.32 

6.55 

2.53 

2.21 

1.15 

> 0.30 

0 

75 

150.15943 

2.29648 

2.34 

0.35 

6.40 

10.40 

1.96 

5.32 

0.22 

1 

76 

150.18268 

2.33601 

2.01 

0.32 

6.31 

-0.32 

1.83 

-0.17 

> 0.55 

0 

77 

150.07148 

2.42307 

3.22 

0.51 

6.28 

3.97 

3.60 

1.10 

> 0.29 

0 

78 

150.09911 

2.40516 

2.51 

0.39 

6.23 

5.50 

2.49 

2.21 

0.46 

0 

79 

150.04249 

2.32799 

2.20 

0.35 

6.19 

3.49 

2.36 

1.48 

> 0.27 

0 

80 

150.13624 

2.26135 

2.06 

0.34 

6.05 

-0.66 

2.16 

-0.30 

> 0.48 

0 

81 

150.12630 

2.41379 

2.21 

0.37 

5.98 

1.27 

2.35 

0.54 

> 0.37 

0 

82 

150.15286 

2.32011 

1.59 

0.27 

5.93 

3.65 

1.60 

2.29 

0.43 

0 

83 

150.02572 

2.31335 

2.75 

0.46 

5.93 

6.50 

2.92 

2.22 

0.42 

0 

84 

150.11186 

2.40879 

2.18 

0.37 

5.92 

-1.58 

2.36 

-0.67 

> 0.46 

0 

85 

150.11984 

2.41767 

2.32 

0.39 

5.87 

7.06 

2.52 

2.80 

0.33 

0 

86 

150.05200 

2.30554 

1.91 

0.33 

5.87 

2.38 

2.23 

1.07 

> 0.28 

0 

87 

150.22409 

2.35646 

3.71 

0.64 

5.83 

-1.10 

3.16 

-0.35 

> 0.59 

0 

88 

150.05389 

2.27630 

2.11 

0.37 

5.68 

0.83 

2.27 

0.37 

> 0.39 

0 

89 

150.16178 

2.26814 

2.15 

0.38 

5.67 

13.81 

2.32 

5.96 

0.16 

1 

90 

150.05476 

2.25801 

2.59 

0.46 

5.63 

-5.87 

3.14 

-1.87 

> 0.41 

0 

91 

150.07011 

2.29022 

1.82 

0.32 

5.60 

3.47 

2.07 

1.67 

> 0.24 

0 

92 

150.05980 

2.40055 

2.37 

0.43 

5.57 

4.56 

3.01 

1.52 

> 0.22 

0 

93 

150.05751 

2.42810 

4.36 

0.78 

5.57 

14.01 

5.47 

2.56 

0.31 

0 

94 

150.06199 

2.37970 

1.95 

0.35 

5.53 

-2.09 

2.40 

-0.87 

> 0.41 

0 

95 

150.01647 

2.32095 

3.42 

0.62 

5.51 

3.28 

3.58 

0.92 

> 0.33 

0 

96 

150.10807 

2.42369 

2.39 

0.45 

5.36 

1.88 

2.82 

0.67 

> 0.32 

0 

97 

150.09548 

2.28661 

1.53 

0.29 

5.31 

0.28 

1.88 

0.15 

> 0.38 

0 

98 

150.16077 

2.34168 

1.54 

0.29 

5.29 

9.34 

1.73 

5.40 

0.17 

1 

99 

150.20984 

2.31258 

2.53 

0.48 

5.29 

15.46 

2.73 

5.67 

0.16 

1 

100 

150.21841 

2.34489 

2.79 

0.53 

5.25 

1.80 

2.72 

0.66 

> 0.38 

0 

101 

150.14854 

2.25458 

2.01 

0.38 

5.22 

-2.33 

2.44 

-0.95 

> 0.41 

0 

102 

150.03720 

2.27215 

2.66 

0.51 

5.21 

14.78 

3.25 

4.55 

0.18 

1 

103 

150.08604 

2.38099 

1.94 

0.36 

5.18 

13.00 

2.28 

5.70 

0.15 

1 

104 

150.14108 

2.42386 

2.29 

0.45 

5.10 

15.70 

2.83 

5.55 

0.15 

1 

105 

150.16471 

2.40932 

2.04 

0.40 

5.04 

4.68 

2.27 

2.06 

0.44 

0 

106 

150.20893 

2.35022 

2.12 

0.42 

5.02 

3.52 

2.40 

1.47 

> 0.25 

0 
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Table A2. The results of the radio/mid-infrared statistical identification process described in Section]^ The columns give the SCUBA-2 
source number, the positions of the adopted optical ID, the VLA 1.4 GHz coordinates (where a radio ID exists), and the relevant mid- 
infrared and radio flux densities, angular offsets (from the SCUBA-2 850 /im position), and corrected probabilities, p, that each association 
could have occurred by chance (given the depth of the supporting data, the relevant number counts, and the counterpart search radius). 
If a given ID is listed more than once, the counterpart with the lowest p-value was treated as a correct association. The robust IDs 
(p < 0.05) are shown in bold, the more tentative IDs (0.05 < p < 0.1) in italics, and the sources for which the optical/near-infrared IDs 
were rejected on the basis of the Zp - zi,\y comparison (see Section |4.2[ | are marked with asterisks. 


ID 

SC 850 - 

R^opt 

deg 

DECopt 

deg 

^^^VLA 

deg 

DECvlA 

deg 

' S ' s.O 

distg.o 

P 8.0 

S24 

mJy 

dist 24 

P 24 

‘ S'VLA 

mJy 

distyLA 

PVLA 

1 * 

150.06460 

2.26405 



18.88 

2.47 

0.062 

0.13 

2.18 

0.036 




2 * 

150.10014 

2.29713 

150.09994 

2.29721 

35.00 

2.23 

0.032 

0.16 

1.42 

0.016 

0.187 

1.85 

0.001 

4 

150.10546 

2.31285 

150.10535 

2.31284 

24.71 

1.57 

0.026 

0.23 

0.65 

0.003 

0.058 

1.62 

0.002 

5 

150.14304 

2.35585 

150.14323 

2.35602 

14.21 

0.66 

0.011 

0.14 

0.40 

0.002 

0.517 

0.20 

0.000 

6 

150.09854 

2.36536 

150.09865 

2.36538 

31.79 

1.35 

0.016 

0.24 

1.20 

0.007 

0.043 

1.60 

0.002 

7 

150.09866 

2.32081 



15.93 

3.16 

0.092 

0.12 

2.30 

0.041 




8 

150.09790 

2.26001 



9.47 

1.91 

0.070 







9 

150.07911 

2.28180 






0.33 

1.35 

0.003 




9 





27.09 

2.21 

0.040 

0.31 

1.97 

0.012 




9 





14.09 

3.85 

0.116 

0.30 

4.05 

0.035 




10 

150.15374 

2.32800 



19.41 

1.36 

0.026 







11 

150.04326 

2.37348 

150.04318 

2.37357 

20.86 

2.15 

0.047 

0.27 

1.30 

0.004 

0.100 

2.03 

0.002 

13 * 

150.08440 

2.29049 



59.24 

2.71 

0.027 

0.42 

2.43 

0.010 




14 

150.10641 

2.25161 

150.10635 

2.25161 

26.43 

2.94 

0.059 

0.58 

2.53 

0.007 

0.112 

2.89 

0.003 

15 

150.11754 

2.32996 



14.97 

1.69 

0.044 

0.16 

1.33 

0.009 




16 

150.05657 

2.37375 

150.05649 

2.37383 

107.95 

0.88 

0.003 

0.46 

0.78 

0.001 

0.088 

0.90 

0.001 

17 

150.20797 

2.38308 



21.93 

0.71 

0.008 

0.07 

0.15 

0.001 




18 

150.16357 

2.37242 

150.16351 

2.37251 

34.47 

1.64 

0.020 

0.56 

1.41 

0.003 

0.138 

1.72 

0.001 

19 

150.11255 

2.37654 



10.47 

0.76 

0.018 

0.10 

0.93 

0.013 




20 

150.15026 

2.36414 



30.45 

1.54 

0.021 

0.21 

1.00 

0.007 




21 

150.09867 

2.31118 



11.97 

0.37 

0.005 

0.18 

0.90 

0.007 




22 

150.05706 

2.29286 



14.36 

2.45 

0.073 

0.16 

2.09 

0.028 




23 

150.12294 

2.36096 



12.93 

0.61 

0.011 

0.23 

0.65 

0.002 




24 

150.10909 

2.29433 



21.57 

1.37 

0.024 

0.22 

0.64 

0.002 




25 

150.03729 

2.34057 

150.03740 

2.34071 

9.31 

2.50 

0.096 

0.07 

1.69 

0.040 

0.062 

1.86 

0.003 

26 

150.07937 

2.34056 

150.07925 

2.34052 

13.55 

3.04 

0.096 

0.16 

2.82 

0.042 

0.061 

3.38 

0.005 

26 








0.18 

2.68 

0.022 




28 

150.12181 

2.34131 



10.69 

1.76 

0.059 

0.10 

2.09 

0.042 




29 * 

150.10525 

2.43499 



13.36 

1.21 

0.030 

0.08 

0.89 

0.014 




31 

150.05248 

2.24555 



37.59 

2.89 

0.043 

0.38 

1.50 

0.003 




33 

150.04098 

2.28063 



11.41 

1.89 

0.062 

0.08 

2.74 

0.068 




34 * 

150.13513 

2.39942 

150.13495 

2.39930 

14.60 

0.38 

0.004 

0.17 

0.32 

0.001 

0.056 

0.96 

0.001 

35 

150.16771 

2.29876 



16.71 

2.85 

0.080 

0.32 

2.62 

0.018 




36 * 

150.08187 

2.41556 






0.64 

0.72 

0.001 




37 

150.06811 

2.27569 



29.14 

1.92 

0.030 

0.45 

1.53 

0.004 




38 

150.07527 

2.37940 






0.13 

1.28 

0.017 




40 * 

150.10531 

2.32590 






0.06 

3.12 

0.050 




42 * 

150.02754 

2.34577 






0.16 

4.32 

0.069 




43 

150.17186 

2.24070 



18.90 

2.87 

0 . 07 A 

0.21 

2.77 

0.032 




44 * 

150.13702 

2.23222 

150.13658 

2.23252 

46.55 

2.98 

0.038 

0.22 

2.39 

0.024 

0.045 

1.94 

0.003 

45 

150.12715 

2.38786 



12.64 

1.21 

0.031 

0.06 

0.98 

0.020 




46 

150.10590 

2.42879 






0.13 

1.99 

0.033 




47 

150.04825 

2.25144 






0.13 

4.21 

0.047 




48 

150.02141 

2.28867 



20.33 

3.16 

0.079 







49 

150.15747 

2.35803 






0.19 

3.38 

0.047 




51 

150.03652 

2.28617 






0.10 

3.63 

0.049 




53 

150.18763 

2.32250 



46.19 

1.80 

0.018 

0.24 

1.73 

0.013 




54 

150.04241 

2.29985 



19.09 

2.46 

0.061 

0.09 

2.47 

0.056 




55 

150.13354 

2.37042 



15.81 

3.32 

0.097 

0.10 

2.10 

0.027 




56 

150.05002 

2.38607 






0.05 

1.15 

0.026 




59 

150.18497 

2.38894 






0.11 

4.48 

0.088 




61 

150.05397 

2.39590 



11.07 

1.39 

0.042 

0.16 

0.73 

0.006 




62 

150.16691 

2.23582 



20.71 

0.89 

0.013 

0.37 

0.69 

0.002 




63 * 

150.07672 

2.39860 



12.09 

2.56 

0.086 







64 

150.13074 

2.31408 






0.18 

4.05 

0.059 




65 * 

150.09156 

2.39904 



58.66 

2.46 

0.024 

0.16 

2.19 

0.018 




66 

150.17561 

2.40159 



27.26 

2.72 

0.052 

0.16 

2.34 

0.034 




67 

150.11132 

2.40320 



9.76 

3.36 

0.121 

0.11 

3.22 

0.067 




68 

150.13001 

2.25269 



10.88 

2.59 

0.092 

0.18 

2.63 

0.034 




71 * 

150.07194 

2.23867 






0.08 

2.14 

0.050 




72 * 

150.06456 

2.32903 



111.50 

2.24 

0.012 

0.29 

1.99 

0.013 




73 

150.20962 

2.35525 

150.20955 

2.35531 

1446.07 

2.32 

0.001 

1.46 

2.07 
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Table A4. The optical-infrared photometry for the SCUBA-2 identifications which lie outside the CANDELS HST imaging. The table gives CFHTLS optical, Subaru ^’-band, 
UltraVISTA near-infrared and IRAC mid-infrared AB magnitudes measured through 3 arcsec-diameter apertures. Errors are the \a values with a minimum adopted error of 0.1 mag. 
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Table A5. The derived physical properties of all 106 sources in the SCUBA-2 sample. The columns show, respectively, SCUBA-2 source 
number, optical spectroscopic redshift (should it exist for a robust ID), optical-infrared photometric redshift Zp, the ‘long-wavelength’ 
photometric redshift zl \ y -, the normalized redshift offset r = (zlw ~ Zp )/(1 + Zp ) (see Section |4.2[ l, a flag indicating the status of the 
redshift information, and our final adopted redshift (with estimated errors), star-formation rate (SFR) and stellar mass (M*). If a source’s 
optical/near-infrared ID was rejected on the basis of an excessive value of r, it is flagged here with 0, if accepted it is flagged with 1 and if 
no ID was found the flag is set to 2. For objects flagged with 1 the final redshift 2 : is the optical-infrared photometric redshift Zp (or Zspec 
if it exists), and therefore a stellar mass can be estimated for the galaxy and is given in the final column. If the flag is 0 or 2, the final 
adopted redshift 2 becomes zl-[y, but no stellar mass can be calculated due to the absence of any optical—infrared photometry. The last 
column gives the source for the spectroscopic redshifts. These include II redshifts from DRI of the zCOSMOS redshift survey undertaken 
in the COSMOS field with the VIMOS spectrograph (Lilly et al. 2007), 2 redshifts obtained as a part of the ST)-HST observations of 
the COSMOS field (Brammer et al. 2012, Skelton et al. 2014) and 2 redshifts from the spectroscopic survey undertaken in the COSMOS 
field with DEIMOS spectrograph (PI: Jehan Kartheltepe). 
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157.8 

± 

10.7 

10.85l°;°i 


29 

0.7270 

n 71 +0.14 
' J^-O.ll 

2 41 +0-24 
^•^1-0.26 

0.97 

0 

2 41 +0-24 
^•^1-0.26 

132.5 

± 

19.5 


VIMOS 

30 


2 (54+0-48 

^■O^-0.41 


2 

2 54+0-48 
^•^^-0.41 

128.9 

± 

11.2 



31 


9 yi'7-|-U.U8 
-0.12 

9 9Q+0.22 

^•^^-0.15 

-0.05 

1 

9 ^^<7+0.08 

^•'l'-0.12 

450.5 

± 

36.5 

11 r)q+u.:^b 

li.ZO_Q 39 


32 



2 02+0-04 
^•'^^-0.45 


2 

2 02+0-94 
^•^^-0.45 

57.8 

± 

7.9 



33 


^•^^-0.65 

3 26"*”®'^^ 
'^•^'^-0.56 

0.25 

1 

2 4n"^^'40 
^•^^-0.65 

206.6 

± 

21.4 

11.66111° 


34 

0.0010 

'^•'^^-0.04 

9 91 +0.41 
^•^1-0.21 

2.21 

0 

9 91 +0.41 
^•^1-0.21 

164.2 

± 

19.1 


3D-ifST 

35 


-*^•'^0-0.16 

^•^ 0 - 0.27 

0.30 

1 

1 3(5+0-24 

l-'^O-0.16 

102.3 

± 

10.5 

10.90lJli 


36 


'^•J^'^-O.ll 

l-'^^-0.14 

1.50 

0 

1 Qn'^0-20 

l-^'J-0.14 

132.5 

± 

19.5 



37 



1 >50+0-20 
l-OO_0.12 


2 

1 56 '^ 0-20 
l-^O-0.12 

261.5 

± 

13.3 



38 


1 08+^-*-^ 
J--yo_0.28 

9 riq+ 0.28 

^•00-0.39 

0.02 

1 

1 Qq + 0.12 

i-yo-0.28 

174.5 

± 

16.8 

10.7ll°;« 


39 



2 (58+0-03 
Z.OO_o.42 


2 

9 cq+ 0.63 
^•^®- 0.42 

356.2 

± 

17.5 



40 


0.8712:27 

^•yy-0.32 

0.65 

0 

^•'J'^-0.32 

447.6 

± 

14.5 



41 



> 3.92 


2 

> 3.92 

171.2 

± 

11.2 



42 

0.9370 

0.96+2'2? 

9 9ri+U.dt) 

^•^y_0.34 

0.65 

0 

9 9ri + U.cJt) 

^•^^-0.34 

299.8 

± 

21.2 


DEIMOS 

43 


ofi+ 0.44 

J--oD_o.26 

9 qq+ 0.79 

^•00-0.49 

0.34 

1 

1 qfi+ 0.44 

l-oO-o 26 

176.2 

± 

22.1 

iu.yu_Q 99 


44 

0.1220 

'^•^'^-0.13 

l.96«;f2 

1.64 

0 

I Q6"^0-38 
l-^°-0.12 

16.8 

± 

2.9 


VIMOS 

45 


q 90+0.47 
‘^•^®-0.13 

3-06lJ:?J 

-0.05 

1 

q 9q + 0.47 
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164.2 

± 

19.1 

10.581J;?® 


46 



1 88+0-21 
-‘-■°^-0.14 


2 

1 qq + 0.21 

4'00-0.14 

132.5 

± 

19.5 



47 


9 t:c+U.V 3 

^•^^-0.45 

9 eq+ 0.63 
^•UO-0 48 

0.02 

1 

9 cc+ 0.75 
^• 0 ^- 0.45 
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± 

32.8 

Il.l2l?;°? 
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Table A5. (continued). 


ID 

^spec 

Zp 

2 LW 

r 

flag 

z 

SFR 

log(M*) 


SC850- 







/Moyr-l 

/M 0 


48 


o 11+0.09 
'^•-^-^-0.16 

^•OO_0.25 

-0.26 

1 

0 11+0.09 
'^•11-0.16 

369.5 ± 37.6 

1 1 qo+0.25 

ii.OZ_o_44 


49 


J--OU_o 20 

^•O0_o.30 

0.17 

1 

4-oy-o.2o 

108.5 ± 11.8 

9.48l°;?; 


50 



l.95t°;«^ 


2 

1 95+0-S7 

J--y'J-0.27 

490.2 ± 15.1 



51 


2 

^•'^^- 0.11 

3.96tl°4 

0.65 

1 

^• 1 ^ 1 - 0.11 

152.4 ±20.8 

10.821®;®® 


52 



o 1 0 + 2.82 

'J-lo_1.28 


2 

0 1 0 + 2.82 
3-l0-i 28 

48.4 ± 8.2 



53 


1 41+014 

1 47+0.09 
-0.08 

0.02 

1 

1 41+0-14 

120.6 ± 9.8 

ll-51+o:S 


54 


o-0'y_0.44 

^•yy-0.30 

-0.24 

1 

'j-yy-0.44 

204.5 ± 22.4 

11 . 021 ?;;® 


55 


1 74+0.11 
-!• '^-0.04 

1 9 C+O .13 

J--^'J_0.08 

-0.18 

1 

1 74+0-11 

J-- '^-0.04 

128.9 ± 11.2 

11 

-‘--‘^•^'-0.16 


56 


9 Qn+0.35 
^•oU_o.4o 

9 1 O+0.48 
^•-^'^-0.31 

-0.18 

1 

9 Qn+0.35 

4 . 0 U -0 40 

194.9 ±23.6 

10.761?;?® 


57 



3 00 + 2.47 
O.OU_o 72 


2 

3 00 + 2.47 

O.OU_Q 72 

143.9 ± 17.0 



58 



2 45 + 0.00 
^•^o_o.49 


2 

2 45 + 0-00 
^•^O_o.4g 

157.8 ± 10.7 



59 


-i-'J^_0.04 

2 04 

^•O^-0.26 

0.13 

1 

1 69"^^'^^ 
J--'J^_0.04 

120.9 ± 13.8 

10.301?;?®5 


60 



> 3.60 


2 

> 3.60 

25.6 ± 3.9 



61 



1 47+0-16 
1-^^ -0.14 


2 

1 4-7+0.16 
1-^'-0.14 

57.8 ± 7.9 



62 


1 ^7+0.08 

2-24t°:i 

0.21 

1 

i.67in? 

129.6 ± 17.3 

10.961°;®® 


63 


n R1 +0.29 

5.4+0;®® 

2.54 

0 

5 . 4 +®;®® 

57.8 ± 7.9 


64 


1 ^0+0.47 

J-.09_o.28 

, t-o-l-0.24 

-0.04 

1 

1 03+0-47 

J--OO-0 28 

85.5 ± 9.4 

10.691?;?; 


65 

2.4750 

Q 1 O+0.07 
O.io_o 03 

z.uy_Q 29 

- 0.11 

1 

2.48 

132.5 ± 19.5 


VIMOS 

66 


-1 yy + 0.08 

4+ '- 0.12 

1 57+0-27 

l-O'- 0.22 

-0.07 

1 

1 177 + 0.08 
!•' '- 0.12 

143.9 ± 17.0 

10.841°'®; 


67 


2 94+0.11 
^•'4^-0.19 

4.yO_Q 23 

- 0.11 

1 

9 04 +O.II 

^•'^^-0.19 

134.3 ± 18.0 

10.091®;®® 


68 


n " 70 + 0.22 
' °-0.18 

i-o'J-o.ie 

0.46 

1 

n 70 + 0.22 
' °-0.18 

16.8 ± 2.9 

9-7ii;;f8 


69 



> 4.59 


2 

> 4.59 

129.6 ± 17.3 



70 



4 07+103 

^•y' -1.79 


2 

4 97+1-00 

^•y' - 1.79 

369.5 ± 37.6 



71 


0 63+012 
U-O9_o 03 

9 ri+O.SO 
^•01-0.44 

1.15 

0 

9 ki+0.80 
^•01-0.44 

450.5 ± 36.5 



72 

2.4460 

u-ou_o 05 

l-OU-0.15 

-0.27 

1 

2.45 

217.9 ± 15.1 

10.40l?'?? 

VIMOS 

73 

0.1660 

n 1 q+0.07 

0 42+0-02 

'4-^^_0.03 

0.22 

1 

0.17 

6.2 ± 0.0 

VIMOS 

74 


2 QQ +010 

^•^^-0.09 

1 94+0-46 

-^•^^-0.31 

-0.26 

1 

2 99+0-10 

^•^^-0.09 

133.2 ± 19.2 

10.891?;;; 


75 


1 "70+0.37 

4-''4_o.13 

9 ri'7+0-23 
^•O'-0.31 

0.12 

1 

1 70 + 0.37 
i-'-^-O-lS 

109.6 ± 13.2 

iU.OZ_o 51 


76 



> 3.58 


2 

> 3.58 

120.6 ± 9.8 



77 


0 64+0 14 
0'-04:_o.09 

2 G1 

^•01-0.65 

1.20 

0 

2 61 + 1-11 
^•Ol_o.65 

57.8 ± 7.9 



78 


9 91 +0.19 

^•^4_o.46 

2 G1 +1-00 
^•01-0.64 

0.12 

1 

9 91 +0.19 
^•^1-0.46 

132.5 ± 19.5 

10.291?;?; 


79 


1 OK+O.OS 
-‘-•^O_o.35 

1 92+0-40 
l-l^^-0.24 

0.30 

1 

1 9[r+0.65 
i-^0_o.35 

48.4 ± 8.2 

Q 00+2.78 

y-oo-i.so 


80 



> 3.42 


2 

> 3.42 

16.8 ± 2.9 



81 


4'-04_o.ii 

> 3.34 

3.35 

0 

> 3.34 

134.3 ± 18.0 



82 


Z.DC)_o.75 


2 

Z.OO_Q 75 

279.8 ± 14.2 



83 


1 9c;+0.05 

4-'40_o.15 

2 04 

^•i'^-0.30 

0.29 

1 

1 ocr + 0.05 

70.5 ± 12.2 

10 31+‘^■22 

iU.Oi_Q 56 


84 

0.3500 

0 3p;+015 

4'-'40_o.io 

> 3.61 

4.19 

0 

> 3.61 

134.3 ± 18.0 


VIMOS 

85 


1 51 +0-26 
1-O1-0.24 


2 

1 K 1 +0.26 

1-O1-0.24 

134.3 ± 18.0 

11 11+0-89 
J--L-J--L_0.52 


86 

1.4530 

1 36+^49 

4-'4o_0.11 

1 qo+ 0-22 

^•'^^-0.15 

-0.03 

1 

1.45 

75.8 ± 10.2 

VIMOS 

87 


2 94+0-01 
^•^^-0.34 

4.82t;;« 

0.48 

1 

2 94+0-01 
^•^^-0.34 

200.9 ±35.4 

10 77+0-85 
J-y-' ' -0.93 


88 


i-'40_o.io 

3.68tlJ® 

1.03 

0 

3.68+;g 

206.6 ±21.4 

Q 70+0.77 


89 

0.9050 

4'-l'0_o.io 

4-DO_q 20 

0.39 

1 

0.91 

25.6 ± 3.9 

3D-HST 

90 





2 

4.yD-o 31 

282.8 ±32.8 



91 


1 41+0.14 

l-^l-0.21 

1 oo+O-O"^ 
J--y'J-0.36 

0.20 

1 

1 41+0-14 

l-^l-0.21 

64.6 ± 11.1 

10 33+0-61 
iU.OO_Q g2 


92 


1 1 9+0.03 

l-l^-0.12 

1 o"7+0-18 

-0.14 

0.12 

1 

1 1 9+0.03 

l-l^-0.12 

57.8 ± 7.9 

io.5il®o;;? 


93 

0.6550 

u.o^_0.09 

3 41 + 107 
'^•^-^-0.85 

1.66 

0 

3 41 + 1-07 
'^•^-^-0.85 

57.8 ± 7.9 


DEIMOS 

94 



> 4.31 


2 

> 4.31 

69.5 ± 5.0 



95 


2.54t°oil 

i.yD_Q 24 

-0.16 

1 

2 54+0-01 
^•O^-0.29 

192.0 ±31.9 

io.67i;;i; 


96 



> 3.12 


2 

> 3.12 

134.3 ± 18.0 



97 



> 3.34 


2 

> 3.34 

171.2 ± 11.2 



98 


0.69l°l® 

1 1O+0.09 

i-i'^-0.08 

0.26 

1 

0.691®;;® 

17.3 ± 2.0 

11 03+1-04 

ii.UCi_Q gi 


99 


u.oy_o.o9 

U.D^_q 04 

-0.07 

1 

0.75 

91.0 ± 2.1 

10.541°;®? 
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> 3.53 


2 

> 3.53 

200.9 ±35.4 



101 



> 3.96 


2 

> 3.96 

16.8 ± 2.9 
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1.7410 

1.0+_o.i4 

1 46+0-^l 

l-^'l-0.16 

- 0.10 

1 

1.74 

133.5 ± 17.5 

10.40l°'j; 

VIMOS 

103 


2 95+0.20 
^•yo_0.55 

1 fi 5 ^+ 0-22 

I-DO-0.28 

-0.32 

1 

2 95+0-20 
^•yo_o.55 

151.8 ±21.6 

10.541?;?? 


104 



1 RR+ 0 - 1 ® 

I.DO-0.29 


2 

1 «q+0-18 

1-O3-0.29 

134.3 ± 18.0 



105 


0 '-O'_ 0.12 

2 09 

^•yy-0.40 

0.85 

0 

^•uy-0.40 

143.9 ± 17.0 
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9 yir+1.66 
^O_o.33 


2 

9 7 [r + 1.66 

^•'0_Q_g3 

6.2 ± 0.0 
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APPENDIX B: MULTI-FREQUENCY IMAGES 
AND SOURCE IDENTIFICATIONS 

In this second appendix we first provide multi-wavelength 
postage-stamp images of all 106 SCUBA-2 sources to illus¬ 
trate the results of the identification process, and then show 
higher-resolution versions of the deepest available near- 
infrared images of the sources, highlighting both the iden¬ 
tified and unidentified sources, along with their estimated 
redshifts. 
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K fl60 SSfim SjjLm 2A(ji.m 450/zm 850/im VLA 


Figure Bl. Multi-wavelength postage-stamp images of all 106 SCUBA-2 sources. Each image is centred on the 850 /i.m source position, 
and is 36 x 36arcsec, with North to the top and East to the left. The blue circle in each stamp shows the search radius adopted for 
each source when searching for counterparts and calculating the significance of any positional association. The yellow and red symbols 
depict the most probable counterparts with p < 0.05 and 0.05 < p < 0.1 respectively. The number given on the 850 ^m stamp is the 
signal-to-noise ratio of the 850 /im detection; the sources are ranked by this number. 
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Figure Bl. (continued). 












































































28 



K /160 S.Gyum 8/xm 24/im 450/zm 850/im VLj4 


Figure Bl. (continued). 




















































































29 



Figure Bl. (continued). 
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Figure Bl. (continued). 
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Figure Bl. (continued) 
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Figure Bl. (continued). 
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Figure Bl. (continued). 
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Figure Bl. (continued) 
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K /160 S.Gyum 8/xm 24/im 450/zm 850/im VLj4 


Figure Bl. (continued). 
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Figure Bl. (continued). 
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Figure B2. Near-infrared (HST Hiqq if available, otherwise UltraVISTA Ks-h&nd) image stamps for all sources in the SCUBA-2 sample. 
Each stamp is centred on the 850 fim position, and is 12 x 12 arcsec, with North to the top and East to the left. The red squares indicate 
the positions of the 450 fim counterparts (when present). The green triangles show the 24 /im ID positions, the yellow circles the 8 ^m ID 
positions, the pink stars the radio ID positions and the red diamonds the final optical ID positions. The colour-coding of the ID names in 
the upper-right corner is as follows. Green means that the adopted optical ID was deemed acceptable on the basis of agreement between 
Zp and zl\y, and the final adopted value of Zp is given in the top-left corner of the stamp. Red means that the proposed optical ID 
was rejected on the basis of inconsistency between Zp and z^wr, and the first number in the upper-left corner of each stamp is then the 
adopted z^y/^ with the second number being the rejected Zp. Where source names are given in blue, this means that no ID was secured 
for the SCUBA-2 source, and only the value of z^y/ was available (and is given in the upper-left corner). Seeveral of the rejected optical 
IDs appear to be rather obvious lenses (SC850-1, 2, 29, 44, 93). 





























































Figure B2. (continued). 
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Figure B2. (continued). 
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Figure B2. (continued). 
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Figure B2. (continued). 








































































Figure B2. (continued). 
































